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1.  PROJECT  SUMMARY 

The  objective  of  this  research  program  is  to  study  theoretically  the  undeilying  principles  of  solid- 
state  dynamics  and  quantum  mechanical  transport  of  carriers  in  ultrasmall  novel  semiconductor  devices. 
The  areas  of  research  are:  1)  theory  of  [rfionon  modes  in  reduced  dimensions,  2)  effects  of  band 
structure  on  electronic  and  optical  properties  of  heterostiucnires,  and  3)  quantum  transport  in  solids 
with  special  emphasis  on  the  non-pertutbadve  role  of  high-electric  fields  and  many-body  effects  in 
dynamical  processes.  The  treatment  of  these  problems  is  mainly  analytical  through  the  development  of 
macroscopic  and  microscopic  physical  models  with  an  emphasis  on  quantum  mechanical  principles.  At 
the  same  time,  numerical  approaches  has  also  bemi  utilized  for  realistic  solutions  with  accuracy. 
Specific  subjects  discussed  in  this  report  include  the  effects  of  confinement  and  localization  on  q>tical 
and  acoustic  {^non  modes,  band  mixing  in  phonon-assisted  transitions,  and  Bloch  electron  quantum 
transport  theory  under  hot-electron  conditions.  The  knowledge  developed  in  this  work  will  be  of  major 
importance  in  explaining  the  novel  {Aienomena  and  fundamental  questions  relating  to  the  breakdown  of 
classical  solid-state  electronics  as  device  dimensional  scales  are  reduced  to  the  submicron  and  ultra- 
sulnnicron  regime. 


2.  PROJECT  DESCRIPTION 


2.1  Introducticm 

Recent  advances  in  semiconductor  materials  giowdi  techniques,  such  as  molecular  beam  epitaxy 
(MBE),  metalorganic  chemical  vapor  deposition  (MOCVD),  and  atomic  layer  epitaxy  (ALE),  have 
made  possiUe  die  fabrication  of  de  '  rith  one  or  more  dimensions  ai^ioaching  the  spacing 
between  planes  of  atoms.  These  nancm'  '-scale  techniques  have  also  opened  new  possibilities  for 
"band  gap  engineering"  of  novel  semiconductor  devices,  including  heteiostnictures  with  spatially 
modulated  energy  band  gaps.  By  using  the  variability  of  ^he  boundary  conditions  which  can  be 
imposed  on  die  wave  functions,  electrical  and  optical  responses  ir  dies?  structures  can  be  tailored  vir¬ 
tually  at  wdL  As  fabrication  technology  has  allowed  such  ultrasmall  strictures  to  be  realized,  many 
new  and  fundamental  questions  have  emerged  concerning  the  underlying  pdiysi(»  of  small  dimensions 
widi  complex,  quantum-scale  boundary  conditions  in  semiconductor  devices.  Important  issues  now 
under  investigation  include  quantum  mechanical  jdienomena  such  as  size  quantization,  {diase  coherence, 
and  highly  ncuiequilibrium  transport  in  which  a  perturbative  treatment  or  the  assurrqition  of  linear 
respmise  is  not  applicable  (see,  for  example.  Refs.  1-3).  More  dian  ever  before,  it  is  important  that  our 
ability  to  analyze  these  physical  phenmnena  occurring  at  ultrasmall  scales  proceeds  technological 
developmoits  and  leads  the  way  to  future  advances. 

In  1990,  the  Office  of  Naval  Research  initiated  sponsorship  of  a  basic  research  program  in  dre 
Department  of  Electrical  and  Computer  En^neering  at  North  Carolina  State  University.  The  general 
goal  of  this  research  program  is  to  study  theoretically  solid-state  dynamics  and  quantum  medianical 
transport  of  carriers  in  ultrasmall  novel  semiconductor  devices.  The  onphasis  is  on  the  developmoit  of 
physical  urKterstanding  of  the  novel  (diermmena  mainly  through  analytical  apfxoaches.  However, 
numerical  evaluation  will  also  be  employed  for  realistic  solutions  with  accuracy.  Our  efforts  have  been 
devoted  to  investigating  three  in^rtant  pioblans.  The  first  one  is  die  dieory  of  phtmon  modes  in 
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heteiDStnictures.  The  changes  in  longitudinal-opdcal  (LO)  phonon  frequencies,  lifetimes,  and  interac- 
tions  witti  carriers  as  a  result  of  reduced  dimensionality  are  the  main  subjects  in  this  part  of  the 
research.  Due  to  the  advent  of  lattice  mismatched  strained-layer  (i.e.,  pseudomoti^c)  stmctures,  the 
effects  of  sti^  on  LO-phonon  modes  have  been  investigated  as  weU.  Recently,  the  scope  of  this 
research  effort  has  been  expanded  to  investigate  the  effects  of  acoustic  phonon  confinement  on 
piezoelectric  scattering  and  deformation  potential  scattering.  The  results  will  be  of  major  importance  to 
a  wide  variety  of  nanoscale  semiconductor  devices  in  which  the  scattering  by  quantized  {^non  modes 
plays  a  significant  and,  at  times,  a  dominant  role  in  determining  the  electronic  and  optoelectronic  pro- 
peities.  The  second  topic  is  the  quantum  mechanical  transport  of  charge  carriers  with  specific 
emi^iasis  on  the  effects  of  electron-electron  and  electron-phonon  interactions  in  the  presence  of  high- 
electric  fields  and  many-body  effects  in  dynamical  processes.  A  novel  formalism  for  treating  Blodi 
electron  dynamics  and  quantum  transport  in  electric  fields  of  arbitrary  strength  and  time  dependence 
has  been  studied  in  an  attempt  to  include  all  quantum  mechanical  effects  collectively  ji  the  lowest 
order  in  the  scattering  strength  [1].  Specific  interest  in  this  inelastic  scattering  problem  arises  from 
issues  and  questitms  relevant  to  the  role  of  high-electric  fields  in  influencing  transport,  electron  relaxa¬ 
tion,  and  noise  generation  as  well  as  ionization  processes  in  quantum  wells,  tunnel  barriers,  superlat¬ 
tices,  atKl  quantum  wires.  At  the  same  time,  attempts  have  been  made  to  study  the  effects  of  quantum 
mechanical  principles  on  carrier  transpoit  with  more  macroscopic  approadies  such  as  dielectric  func¬ 
tion  theory.  As  the  third  topic,  we  have  recently  begun  to  emphasize  the  effects  of  band  structure 
which  are  critical  to  understanding  electrortic  and  optoelectronic  properties  of  heterostructures. 
Detailed  knowledge  of  band  structure  in  the  presence  of  heterointerfaces,  sudi  as  mixing  and  evanes¬ 
cent  states,  is  essential  for  the  characteristic  response  of  a  mesoscopic  system  whidi  is  beyond  the 
realm  of  the  conventional  effective  mass  {qrproximation.  The  main  subjects  of  interest  in  this  topic 
include  tiie  effects  of  band  mixing  in  linear  and  nonlinear  optical  constants,  timneling,  and  phormn- 


assisted  transitions. 
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Although  only  four  years  old.  this  program  has  already  resulted  in  thiity  nine  refereed  publica¬ 
tions  in  the  literature,  five  additional  manuscripts  are  currently  in  press,  and  one  more  has  been  submit¬ 
ted  to  a  major  technical  journal.  Numerous  invited  talks  and  presentations  have  been  and  will  be  given 
at  conferences  and  woiicshops  throughout  the  United  States  and  in  other  countries.  A  listing  of  these 
publications  is  given  in  Appendix  A.  In  addition,  this  research  program  provides  a  perfect  complement 
to  other  aspects  of  our  research  efforts  in  which  numerical  modeling  aspects  of  transport  study  are 
emphasized.  Due  to  the  complicated  geometry  and  boundary  conditions,  a  realistic  model  for  a  realis¬ 
tic  device  structure  is  attainable  only  through  numerical  approaches.  The  accurate  i^ysical  models  and 
understaiKling  developed  in  this  research  facilitates  the  development  of  tiie  numerical  models  better 
able  to  explain  experimental  observations  and  to  predict  new  physical  phenomena.  Overall,  this 
research  program  supported  by  ONR  has  been  efficiofit  and  productive.  The  quality  of  our  program 
will  continue  to  improve  in  the  future.  In  this  annual  report,  the  progress  and  accomplidiments  made 
during  die  past  contract  period  will  be  summarized. 

The  structure  of  the  rest  of  this  report  is  as  follows:  in  Sectitm  2.2,  the  researdi  results  of  the 
past  contract  period  are  summarized;  Section  2.3  provides  the  list  of  publications  resulting  from  this 
research  supported  by  ONR  during  the  1993-94  contract  period;  Section  3  contains  background  infor¬ 
mation  describing  project  personnel;  and  the  Aj^ndices  include  a  list  of  refereed  publications  sup¬ 
ported  by  this  research  program  since  1990  and  the  tide  page  of  each  paper  published  during  the 
1993-94  contract  period. 
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2.2  Summary  of  Research  Results 

This  section  provides  the  current  status  of  our  research  efforts  on  the  theory  of  electron-acoustic- 
{^non  interactions,  the  LO  phonon  lifetime,  the  effects  of  band  mixing,  and  quantum  mecharucal  tran¬ 
sport  in  semiconductor  nanostructures.  The  major  accomplishments  during  the  past  contract  period  are: 

•  We  have  derived  ai^priately  normalized  tq)proximate  expressions  for  die  acoustic  [dionons 
confined  in  a  free-standing  rectangular  quantum  wire  by  quantizing  the  acoustic  {dionon  displace¬ 
ments  and,  subsequently,  have  formulated  the  interaction  Hamiltonian  for  the  deformation  potoi- 
tial  associated  with  confined  acoustic  phonon  modes  in  rectangular  quantum  wires.  The  defoima- 
timi  potential  scattering  rates  due  to  the  confined  acoustic  phonon  modes  have  been  obtained  for 
GaAs  quantum  wires  with  a  range  of  cross  sectional  dimensions.  Comparison  with  the 
corresponding  rates  obtained  with  bulk  phonons  demonstrates  diat  a  proper  treatment  of  confined 
acoustic  {dionons  may  be  essential  to  correctly  modeling  electron  scattering  rates  at  low  energies 
in  nanoscale  structures. 

•  We  have  develt^ied  an  approach  to  calculate  the  lifetime  of  LO  phonons  via  emission  of  two 
acoustic  phonons  in  bulk  zincblende  semiconductors.  The  interaction  Hamiltonian  for  anhar- 
monic  decay  is  derived  based  on  Keating’s  treatment  of  anharmonic  contributions  in  the  elastic 
strain  oiergy  of  a  crystal  Application  of  this  model  to  bulk  GaAs  provides  excellent  agreement 
with  available  experimoital  data.  Since  the  parameters  employed  in  the  model  can  be  obtained 
erqwrimentidly  and  theoretically,  our  aiq)roach  provides  a  useful  tool  for  investigating  LO  [tenon 
lifetimes  in  semicmiductors. 

•  We  have  calculated  the  F-X  intersubband  scattering  rates  due  to  optical  [tenon  emission  by 
using  a  tight  Unding  method  for  electrons  and  the  dielectric  continuum  model  for  optical  [te¬ 
nons,  and  have  demonstrated  that  this  is  an  important  mechanism  for  the  F-X  relaxation  process 
in  type-n  superlattices.  Among  the  various  phonon  modes,  die  strongest  are  found  to  be  die 
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AlAs  confined  modes  and  the  weakest  are  the  GaAs  confined  modes.  It  is  also  shown  that  the 
parity  of  election  states  affects  the  transition  rates  significantly.  The  theoretical  results  olHained 
from  our  calculations  are  in  general  agreement  with  the  experimental  data. 

•  We  have  developed  an  extension  of  the  Knxiig-Penney  model  to  treat  mixings  of  different  valleys 
in  type-II  superiattices.  The  value  of  the  mixing  parameter  used  in  this  model  was  deduced  by 
comparison  with  more  sophisticated  theory  and  various  experiments  for  dw  GaAs/AlAs  material 
system.  It  is  shown  that  with  the  use  of  a  single  value  for  the  mixing  parameter,  this  simple  and 
efficient  model  can  reasonably  reproduce  the  important  aspects  of  valley  mixing  effects  in 
GaAs/AlAs  superiattices  over  a  wide  range  of  layer  thicknesses. 

•  We  have  developed  a  novel  fonnalism  for  treating  Bloch  electron  dynamics  and  quantum  tran¬ 
sport  in  a  superimposed  uniform  and  oscillatory  electric  field.  Based  on  this  treatment,  a  com¬ 
plete  analysis  of  single-band  and  multi-band  {xxxxsses  and  the  total  transition  rates  including  all 
possible  transitions  have  been  obtained.  We  have  also  investigated  the  effea  of  localized  impuri¬ 
ties  on  transport  of  an  electron  in  a  unifonn  electric  fidd,  and  have  formulated  a  fiieoiy  to 
describe  Blodi  electron  transport  through  spatially  localized  inhomogeneous  potential  barriers 
using  a  field-dependent  Green’s  function  approach. 

Each  of  diese  results  is  presented  in  detail  in  the  following  discussions. 

Electron-confined-acoustic-rtKMion  scattering  in  quantum  wires 

A  number  of  proposed  plications  of  mesoscopic  electronic  structures  involve  carrier  transport  at 
low  tmnperatures  and  low  carrier  energies;  frequently,  the  regime  of  interest  is  one  where  dimoisional 
confinement  modifies  the  phase  space  substantially.  It  is  well  known  that  in  this  low  tmnperature,  low 
energy  regime,  acoustic  phonrms  play  an  enhanced  role  in  carrier  scattering  and  may  dominate  over  the 
scattering  of  carriers  by  optical  i^nons.  In  addition,  in  nanoscale  structures  it  is  possible  that  phase 
space  lestrictitHis  may  weaken  or  forbid  optical  i^KHron  scattering  processes  that  would  normally 
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dominate  in  bulk  stnictures.  Recently,  there  has  been  an  extensive  effort  to  study  the  role  of  dimen¬ 
sional  confinement  in  modifying  LO  phonon  modes  and  their  interactions  with  charge  carriers  in 
nanoscale  and  mesoscopic  semiconductor  structures,  as  demonstrated  by  numerous  publications  by  us 
(supported  by  this  program)  as  well  as  by  others.  On  the  other  hand,  there  are  relatively  few  treat¬ 
ments  dealing  with  the  role  of  dimensional  confinement  in  modifying  acoustic  phonon  modes  and  their 
interactions  with  charge  carriers.  In  spite  of  the  fact  that  there  is  an  extensive  literature  on  the  theory 
of  acoustic  modes  in  conventional  waveguides,  resonators  and  related  structures,  few  efforts  have  been 
reported  to  formulate  a  theory  of  acoustic  phonons  in  nanoscale  structures  where  both  phonon 
confinement  and  a  quantum  mechanical  treatment  of  i^onon  normalization  are  essential.  The  necessity 
for  such  theoretical  treatments  has  been  demonstrated  recently  by  experimental  studies  providing  both 
direct  and  indirea  evidence  [4,S]  of  the  importance  of  acoustic  phonon  confinement  in  reduced  dimen¬ 
sional  electronic  structures. 

During  the  past  contract  period,  we  have  initiated  a  research  effort  to  study  the  effects  of 
confinement  on  acoustic  phonons  and  their  interactions  with  electrons  in  nanoscale  structures.  As  a 
first  project  in  this  area,  we  have  developed  appropriately  normalized  approximate  expressions  for  the 
acoustic  [^nons  confined  in  a  free-standing  rectangular  quantum  wire  by  quantizing  die  acoustic  [dio- 
non  displacements  and,  subsequently,  have  formulated  the  interaction  Hamiltonian  for  the  deformation 
potential  associated  with  confined  acoustic  [dionon  modes  in  rectangular  quantum  wires.  In  addition, 
preliminary  works  have  been  done  to  determine  the  role  of  acoustic  phonon  confinement  in  modifying 
piezoelectric  scattering  as  well  as  deformation  potential  scattering  in  a  number  of  different  quantum- 
wire  and  quantum-dot  geometries  of  significant  interest  in  nanoelectrcmics  and  optoelectronics.  In  this 
section,  we  present  our  study  on  die  effects  of  confinemem  on  acoustic  {dionon  modes  in  a  free¬ 
standing  rectangular  quantum  wire  and  the  associated  deformation  potential  scattering. 
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The  specific  stracture  of  our  interest  is  a  finee-standing  rectangular  quantum  wire  of  infinite  length 
in  die  z-directitm  having  an  x-directed  height  (or  thickness),  2a.  and  a  y-directed  width.  2d;  the  origin 
of  coordinates  in  the  x-y  plane  is  taken  to  be  at  the  geometric  center  of  the  rectangular  cross  section. 
Our  appioadi  on  confined  acoustic  phonon  modes  is  based  on  the  classic  study  by  Morse  on  acoustic 
modes  in  a  convendtmal  waveguide  of  rectangular  cross  section  [6].  Through  an  extended  analysis. 
Morse  has  derived  an  a{^roximate  set  of  hybrid  compressional.  or  dilatational,  acoustic  modes  which 
are  found  to  accurately  approximate  the  experimentally  observed  modes  over  a  wide  range  of  condi¬ 
tions.  Specifically,  he  has  found  that  the  i^roximate  hybrid  modes  derived  by  assuming  separate 
boundary  ctmditions  have  simple  analytical  representations  (as  in  our  previous  analysis  of  optical  i^io- 
non  modes  in  rectangular  quantum  wires)  and  provide  convenient  a^roximations  for  the  rectangular 
geometry  when  the  cross-sectional  dimensions  have  aspect  ratios  of  approximately  2  or  greater. 
Adopting  Morse’s  form  for  the  approximate  separation-of-variable  solution,  the  x-,  y-,  and  z-directed 
acoustic  mode  displacements  for  the  "thickness"  mode  may  be  written,  respectively,  as 

Ui  =  A{sink]X  +  asink2x}cos(hy)e^*^^ ,  (la) 

vj  =  A{-^coskiX  +  pcosk2x}sin(hy)e^*^) ,  (15) 

wi  =  jA{— ^coskjx  +  -^-(k2a  +  hp)cosk2x}cos(hy)e^*“**^ ,  (Ic) 

*1  T 

where  A  is  the  normalization  constant,  y  is  the  z-directed  free  wavevector,  c  is  die  phase  velocity,  a 
and  P  are  defined  by 


sinkia  2(h^  +  y^) 
sink2a  (V^  +  h^-  k^)  ’ 


sinkia  2k2h 
sink2a  (-y2  +  52  _  ^2) 


(2a) 

(2b) 


and 


k?  +  h2  =  7*Uc/Cd)2-ll. 

k|  +  h*  =  72Uc/cj2-ii. 
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(3a) 
(3b) 

with  the  compressional.  or  dilatational.  sound  speed.  Cj.  and  the  transverse,  or  shear,  sound  speed,  c,. 
The  boundary  oonditi<Mis  for  acoustic  inodes  in  hee-standing  structures  require  that  the  stress  com¬ 
ponents  normal  to  die  surface  vanish  at  the  surface.  Application  of  this  condition  at  xs±a  (i.e.. 
^  r^^)  results  in  an  equation  which  serves  as  the  dispersion  relation. 

tank2a  4kik2(h^  +  Y*) 
tank^a  "  (h^  +  7*  - 

However,  the  corresponding  boundary  condition  at  y=±d  (i.e..  Ty^T^ysT^^)  cannot  be  satisfied 
exactly  due  to  the  well-known  difficulties  at  the  comers.  When  d>2a.  this  problem  can  be  circum- 
voued  since  the  two  shear  stresses  T;(y  and  become  negligiUe.  Accordingly,  only  the  extensional 
stress  Tyy  needs  to  vanish  at  the  surface:  this  condition  requires, 

h  =  (n  +  .  n  s=  0,1.2,  •  •  •  .  (5) 

2  d 

The  imncipal  thickness  mode  (i.e.,  tt=0  or  h=ic^)  has  no  nodal  surfaces  parallel  to  the  length,  and  is 
the  dominant  mode.  For  a  given  h  and  y,  both  k|  and  k2  may  be  written  in  terms  of  c  as  shown  in  Eq. 
(3).  When  substituted  in  Eq.  (4),  the  jdiase  velocity  c  has  multiple  solutions  due  to  the  periodic  nature 
of  the  trigonometric  functions.  Cmisequently,  the  dispersion  relatitm  or  die  phonon  finequency  axy  (=cy) 
in  a  quantum  wire  needs  an  additiorul  index  m  to  distinguish  differem  modes. 

In  addition  to  the  thidmess  modes,  another  set  of  acoustic  modes  is  observed  experimentally  [6]. 
These  modes  corresponds  to  "widdi  modes"  and  are  determined  in  a  manner  similar  to  that  used  to 
ddermine  die  diickness  modes.  By  satisfying  the  boundary  conditimis  (»i  the  stress  at  y=±d  (i.e.. 
Ty^T^T^xO),  die  solutions  for  the  widdi  modes  show  expressions  analogous  to  Eq.  (1)  widi  die 
roles  of  X  and  y  as  well  as  k  and  h  interchanged,  respectively.  For  diese  modes,  k  is  then  determined 
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by  approximate  boundary  conditions  at  x=±a.  The  dispersion  relation  for  the  width  mode  is  identical 
in  form  to  Eq.  (4).  As  for  the  thickness  modes,  the  principal  mode  widi  k=0  is  die  dominant  mode. 
Figure  1  deincts  the  dispersion  curves  of  the  six  lowest  thickness  modes  (h=ic/2d,  m=l,  ■  •  ■  ,6)  and  the 
corresponding  width  modes  (k=0,  m=:l.  ■  •  ■  .6)  for  a  GaAs  quantum  wire  free-standing  in  vacuum. 
The  quantum  wire  cross  sectional  dimension  is  chosen  to  be  28.3  A  by  56.6  A  (i.e..  10  by  20  mono- 
layers).  As  expected,  the  width  modes  tend  to  have  lower  energies  than  the  thickness  modes  since  the 
width  is  greater  than  the  thickness  for  the  case  represented  in  this  figure. 

When  properly  quantized,  the  acoustic  modes  derived  in  conventional  waveguides  can  be  applied 
to  describe  the  characteristics  of  confined  acoustic  phonon  modes  in  semiconductor  quantum  wires. 
Following  the  quantization  procedure,  the  nonnalization  constant  in  Eq.  (1)  may  be  determined  by 

^/dx  /dy  {u,ur  +  V,,;  +  w,w,-}  -  .  (6) 

where  M  is  the  mass  for  a  unit  cell  and  tOy  is  the  finequency  corresponding  to  the  mode  with  wavevec- 

tor  Y  in  the  z  direction,  as  described  before.  With  this  condition,  the  classical  phorxm  displacemoit 
may  be  written  in  terms  of  quantized  normal  coordinates,  and  the  interaction  Hamiltonians  caused  by 
lattice  vibrations  can  be  developed  readily.  In  particular,  the  interaction  Hamiltonian  for  deformation 
potential  scattering  is  given  by, 

'  « 

Hdrf  =  E.V.TO)  =  E.2  [c„.™(Y)  +  <^(-^)]  ^  +  ^  +  ^  .  (7) 

where  c^^Cy)  and  are  the  usual  annihilation  and  creation  operators.  In  this  equation,  die  sum 

over  Y  represents  the  usual  integration  over  wavevector,  while  the  sums  over  n  and  m  represent  the 
addition  of  the  various  acoustic  phonon  modes.  When  the  expressions  given  in  Eqs.  (1)  and  (6)  are 
adopted  for  determining  the  phonon  displac^ents,  this  Hamiltonian  provides  die  desired  descri{ttion 
for  the  electron-acoustic-phonon  interaction  in  rectangular  quantum  wires,  and  replaces  the  commonly- 
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Hgure  1.  Dispersion  curves  for  the  six  lowest  width  and  thickness  modes  (m  =  1,  *  *  ■  .6)  of  a  28.3 
A  X  56.6  A  GaAs  quantum  wire.  The  solid  lines  are  for  the  width  modes,  and  the  dashed  lines  are  for 
the  thickness  modes. 
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Figure  2.  Defonnation  potential  scattering  rales  for  bulk  and  confined  acoustic  [^non  modes  in  a 
28.3  A  X  56.6  A  GaAs  quantum  wire  at  77  K.  Enhancements  in  the  scattering  rates  for  die  case  of 
ctmfined  acoustic  modes  occur  at  the  onset  of  emission  for  the  various  width  and  thidcness  modes. 
These  diresholds  are  at  0.03,  2.36,  2.55,  4.90,  7.30,  and  7.40  meV  for  die  width  inodes,  and  at  2.06, 
4.44,  5.90,  9.87,  14.5,  and  15.1  meV  for  the  diickness  modes.  The  {dotting  resolution  dqricted  is  not 
fine  enough  to  illustrate  folly  the  im{)ortance  of  the  density-of-states  effects  in  the  rectangular  quantum 
wire. 
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used  Hamilttmian  based  on  bulk  acoustic  i^nons. 

Based  on  die  deformation  potential  interaction  Hamiltonian  derived  above,  we  have  calculated  die 
electron-acoustic-plKinon  scattering  rates  in  a  GaAs  rectangular  quantum  wire  free-standing  in  vacuum. 
For  simplicity,  an  isotropic  cubic  medium  has  been  assumed  and  die  compressional,  or  longitudinal, 
sound  speed  has  been  taken  to  be  that  of  GaAs;  it  should  be  noted  that  imposing  both  of  these  con¬ 
straints  makes  it  impossible  to  have  a  transverse  sound  speed  matching  that  of  GaAs.  This  is  a  conse¬ 
quence  of  the  fact  diat  GaAs  may  be  treated  as  having  an  isotropc  elastic  tension  only  in  a  very  rough 
aiqiroximatioiL  In  this  work,  Poisson’s  ratio,  o,  is  taken  to  have  the  {diysically  reasonable  value  of 
1/3;  diis  choice  fixes  the  value  of  c,  for  a  given  c^  in  Eq.  (3).  The  Fermi  Golden  Rule  scattering  rates 
for  emission  and  absorption  are  plotted  in  Fig.  2  as  funr^ons  of  electron  energy.  In  diis  figure,  die 
extreme  quantum  limit  is  assumed  and  the  rates  are  obtained  for  the  ground  (electronic)  states  at  77  K. 
For  comparison,  the  results  for  bulk  acoustic  modes  are  shown  as  welL  Two  distinct  and  important 
features  are  obvious  from  Fig.  2.  Hrst,  the  appearance  of  stracture  is  prominent  in  die  scattering  rates 
(for  confined  phonons)  which  results  from  the  energy  threshold  for  the  different  mode  values,  m,  of  the 
thickness  and  width  modes.  As  can  be  seen  from  Fig.  2,  each  of  these  modes  makes  a  notable  contri¬ 
bution  to  the  density-of-states  and  to  the  scattering  rate.  In  particular,  the  scattering  rates  at  tow  ener¬ 
gies  show  pronounced  peaks  and  are  strongly  enhanced  due  to  the  dominance  of  selected  compres¬ 
sional  modes  in  die  emission  process.  Thus,  it  is  essential  to  retain  a  number  of  acoustic  modes  for  an 
accurate  estimation  of  scattering  rates.  The  results  shown  in  Fig.  2  have  beoi  obtained  by  including 
the  six  lowest-older  thidoiess  modes  as  well  as  the  six  lowest-order  width  modes.  Due  to  die  limited 
lesolutitMi  in  {dotting,  die  details  of  the  one  dimensional  nature  (such  as  die  number  of  {leaks  and  dieir 
heights)  are  not  illustrated  fully  in  this  figure.  The  second  im{)oitant  feature  of  Fig.  2  is  that  the 
scattering  rates  for  die  case  of  the  confined  acoustic  modes  are  higher  than  the  cone^nding  scattering 
rates  based  on  bulk  {dunons.  These  enhanced  scattering  rates  {irovide  an  indication  that  conceptual 
designs  for  mesosco{dc  devices  need  to  be  based  on  an  awareness  of  the  faa  diat  confined  acoustic 
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modes  may  play  a  significant  role  in  carrier  transport  in  these  devices.  It  is  also  important  to  note  that 
the  sets  of  the  solutions  for  confined  acoustic  modes  given  above  (i.e.,  thickness  and  width  modes) 
may  imt  be  comi^ete.  Thus  it  is  plausible  that  confinement  may  enhance  the  electron-acoustic-phonon 
scattering  rate  even  more  than  it  is  presented  in  this  work. 

Theoretical  calculation  of  LO  phonon  lifetimes 

In  polar  semiconductors,  carrier  energy  and  momentum  relaxation  processes  are  frequently  dcxn- 
iruited  by  the  interaction  of  carriers  widi  optical  phonons,  particulariy  the  LO  phonmis.  It  is  now 
well-known  that  carrier  relaxation  rates  underlie  a  range  of  phenomena  other  than  simple  carrier  tran¬ 
sport  Previous  analyses  of  carrier  relaxation  in  semiconductor  microstiuctuies  have  demonstrated  that 
emission  alone  does  not,  in  general,  dominate  the  carrier  energy  loss  rates  since  subsequent  (^cal 
idionon  absorption  events  slow  the  relaxation  process.  Instead,  it  has  been  shown  that  the  overall  car¬ 
rier  relaxation  rate  is  frequently  domiiuited  by  the  decay  of  the  strongly  interacting  optical  phonon  into 
weakly  interacting  acoustic  phonons.  Thus,  the  detailed  knowledge  of  optical  {rihontHi  decay,  which  is 
characterized  by  lifetime,  is  of  major  importance  in  understanding  carrier  dynamics  in  soniconductors. 
The  decay  process  of  t^tical  i^nons  arises  primarily  from  die  diree-idionon  interaction  throu^  the 
anharmonic  terms  of  the  crystal  potential  energy.  This  phenomemm  has  been  scrutinized  by  a  number 
of  authors  mostly  in  bulk  materials  [7-12].  The  approaches  taken  by  these  audwrs  can  be  classified 
broadly  in  two  categories:  highly  complex  and  complete  microscopic  models  [7,8,11]  with  parameters 
which  are  very  difficult  to  measure,  and  simpler  macroscopic  treatments  [9,10]  where  attonpts  have 
beat  made  to  approximate  various  anharmonic  contributions  by  replacing  them  with  a  single  parameter 
related  to  an  average  of  third  order  elastic  constants  or  a  Griineisoi  constant  An  approach  similar  to 
the  latter  treatments  has  been  adopted  recently  to  estimate  optical  [diooon  lifetimes  in  quantum  wdls  by 
taking  into  account  the  oivelope  functions  of  the  confined  optical  {dionons  [12].  A  major  difficulty 
associated  with  die  simple  approaches  is  that  drey  rely  on  finding  a  valid  eiqiression  for  this  appropri- 
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atefy  averaged  single  parameter  from  the  experimentally  measuraUe  quantities.  Otherwise,  it  simply 
acts  as  a  fittii^  parameter  without  a  solid  physical  basis. 

Last  year,  we  began  an  investigation  on  anharmonic  decay  of  LO  phonons  in  nanostnictures.  The 
principal  aim  of  this  effort  is  to  study  the  effects  of  reduced  dimensions  on  the  LO  i^non  lifetime 
and  the  resulting  carrier  relaxation  processes  throu^  the  development  of  a  {riiysically  valid  model  To 
facilitate  the  capability  to  predict,  the  sought-after  model  must  incorporate  parameters  based  on  a  firm 
I^ysical  footing,  which  can  be  either  measured  experimentally  or  estimated  from  odrer  theoretical  ana¬ 
lyses.  Sudi  a  model  without  fitting  parameters  has  not  yet  been  formulated  even  in  bulk  materials  not 
to  mention  in  heterostmctures.  As  a  first  step  toward  this  goal,  we  have  developed  an  approach  to  cal¬ 
culate  the  lifetime  of  LO  phonons  via  emission  of  two  acoustic  phonons  in  bulk  zincblende  sonicon- 
ductors.  The  interaction  Hamiltonian  for  anharmonic  decay  is  derived  based  on  Keating's  treatmem  of 
anharmonic  contributions  in  the  elastic  strain  energy  of  a  crystal  [13].  Application  of  diis  model  to 
bulk  GaAs  provides  excellent  agreement  widi  available  experimental  data.  Since  the  parameters 
employed  in  the  model  can  be  obtained  experimoitally  and  theoretically,  our  tqrproadi  provides  a  use¬ 
ful  tool  to  investigate  LO  phonon  lifetimes  in  semiconductors.  A  brief  summary  of  the  formalism  and 
results  of  our  calculations  are  provided  below. 

Whnr  a  crystal  potential  is  expanded  in  powers  of  displacements  of  the  atoms  from  dieir  equili¬ 
brium  positions,  we  obtain  a  quadratic  term  along  with  cubic,  quarlic  and  other  higher  order  terms. 
For  simple  analysis  of  foe  dispersion  relation,  it  is  acceptable  to  consider  only  the  quadratic  term  and 
ignore  the  rest  of  foe  higher  order  terms  in  what  is  known  as  the  harmonic  qrproximation.  However, 
the  harmonic  tqrproach  cannot  describe  the  decay  of  phonon  modes  caused  foe  cuUc  and  other 
higher-order  terms,  also  known  as  the  arfoarmonic  terms,  in  the  crystal  potential  It  is  generally 
accepted  that  tiie  cubic  term  dominates  over  all  arfoarmonic  terms  in  ffoonon  decay.  For  this  process 
involving  three  ffoonons,  the  interaction  Hamiltonian  may  be  written  as 
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where  ic,  1^,  and  1^'  0>  j^  ^  j")  represent  the  phonon  wavevectors  (polarization  modes).  N  is  the 
number  of  unit  cells,  and  P  describes  the  cubic  coupling.  The  displacement  for  the  phonon  mode  in 
normal  coordinates  may  be  represented  as 

f  11/2  ^  ^ 

where  aj^  and  aj^  are  the  annihilation  and  creation  operators,  is  die  polarization  vector,  m  is  the 
average  mass  of  the  lattice  atoms,  and  is  the  hequency  of  the  normal  mode.  For  LO  [dionon 
modes,  die  decay  occurs  mainly  through  the  creation  of  two  longitudinal-acoustic  (LA)  modes.  At  die 
same  time,  die  normal  process  dominates  over  the  umklapp  process  whoi  the  LO  phoncm  wavevector 
is  srruill.  Retaining  only  this  term,  the  decay  Hamiltonian  for  LO  phonons  becomes 


The  indices  for  mode  polarization  have  been  removed  since  only  (me  decay  path  has  been  considered 
(i.e..]cix,->lc'i^  +  1^'la)-  Contributions  by  other  processes  can  be  formulated  easily  hom  Eq.  (8). 

In  the  interaction  Hamiltonian  derived  above,  all  the  terms  odier  than  are  known. 

However,  it  is  tpiite  difficult,  if  not  impossible,  U)  properly  measure  die  anharmonic  potential  term.  To 
circumvett  this  problem,  Keating  [13]  has  used  the  theory  of  elasticity  since  diere  exists  a  relationship 
between  the  dtird  order  elasticity  coefficients  and  the  anharmonic  term  in  the  crystal  potential.  Elasti¬ 
city  coefficients  also  are  related  closely  to  tire  Griineisen  parameter.  From  the  analysis,  the  strain 
energy  densiQr  U,  associated  with  the  third  order  contributitms  in  cultic  crystals  reduces  to  [13] 

U,  =s  -^Cin(ei  +  el  +  63)  +  ■|■Cll2[e|(e2  +  C3)  +  e2(e3  +  ej)  +  t^(ci  +  ej)]  +  C]23eie2e3  (ii) 
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+  yCi44(®ie4  +  ®2^  +  C3C6)  +  +  Cg)  +  t^el  +  C4)  +  63(64  +  65)] 

+  C4S6e4e5®6  +  -^^11(6?  +  62  +  63)  +  -jCi2[6?(62  +  63)  +  62(63  +  61)  +  63(61  +  62)]  . 

wh6ie  Cjkl  (Cjk)  are  th6  third  (s6Cond)  ord6r  elastic  constants,  respectively,  and  ej  represents  the 
linear  part  of  die  strain  variable  as  defined  in  Ref.  13.  A  careful  comparison  of  Keating’s  fonnulation 
with  that  of  a  microscopic  model  by  Tua  and  Mahan  [11]  reveals  that  U,  corresponds  to  P  as  defined 
in  Eq.  (8);  accordingly,  a  valid  expression  for  the  decay  Hamiltonian  can  be  developed  once  U,  is 
evaluated,  hi  general,  the  values  for  the  second  order  and  third  order  elastic  constants  can  be  obtained 
from  the  literature.  At  the  same  time,  the  magnitude  of  the  strain  variable  ej  may  be  estimated  based 
on  numerical  models  for  lattice  dynamics  such  as  the  valence  shell  model.  Based  on  the  Hamiltmiian 
derived  above,  the  LO  phonon  decay  rate  can  be  calculated  readily  by  following  Fermi’s  golden  rule 
appioximatioiL  Then,  the  lifetime  for  the  LO  phonon  witii  wavevector  ^  may  be  defined  as 

where  n;*  is  the  normal  Bose-Bnstein  phonm  occtg)ati(xi  number.  When  more  than  one  path  exists  for 
decay  (i.e.,  patiis  other  than  L0->LA4LA).  eadi  contribution  is  added  to  obtain  the  LO  j^ioncHi  life¬ 
time.  In  this  case,  the  sum  in  Eq.  (12)  is  over  both  the  phonon  wavevector  and  mode  polarization. 

Using  this  formalism,  we  have  calculated  tire  lifetimes  of  bulk  LO  i^ionons  in  GaAs.  In  these 
calculations,  tiie  values  for  ej  have  bear  estimated  based  on  our  previous  work  on  tire  valoK»  shdl 
model  (whidi  was  discussed  in  detail  in  our  rqport  last  year).  Although  elastic  omstants  fw  GaAs  are 
available  in  the  literature,  stnne  of  die  measured  third-order  constants  (C123  and  €144  in  particular) 
show  wide  variations.  To  circumvent  this  uncertainty,  we  have  used  a  rdationship  between  die  third 
order  elastic  omstants  and  die  microscopic  force  constants  developed  by  Keating  [13],  and  have  diosoi 
an  optimum  set  systnnatically.  Ihe  values  adopted  in  diis  study  are  listed  in  Table  1.  These  values 


Lifetime  (10 


Table  I.  Tbiid-oider  elastic  constants  for  GaAs  in  the  unit  of  GPa. 


Figure  3.  LO  i^non  lifetime  as  a  function  of  wavevector  in  bulk  GaAs  along  the  <100>  and  <111> 
directions  at  77  K  and  300  K. 
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are  well  within  the  range  of  experimoital  variations.  Figure  3  shows  the  LO  {rfionon  lifetime  as  a 
function  of  wavevector  in  bulk  GaAs  along  the  <100>  and  <1 1 1>  directions  at  77  K  and  300  K.  Since 
the  decay  into  two  LA  phonons  is  the  dominant  process,  we  have  considered  only  this  mechanism. 
The  other  processes  are  either  not  possible  (when  energy  and  momentum  cannot  be  conserved)  or  pro¬ 
vide  a  small  contribution  to  the  overall  lifetime.  As  can  be  seen,  the  lifetime  decreases  widi  increasing 
temperature  maiidy  due  to  the  Bose-Einstein  occupation  number.  The  results  obtained  horn  our 
analysis  agree  well  with  the  widely  accepted  experimental  data  (3.5  ps  at  300  K  and  7  ps  at  77  K). 
The  agreement  is  better  at  ro(»n  temperature  than  at  low  temperatures.  This  is  reasonable  since  the 
third-order  elastic  constants  used  in  our  study  are  estimated  at  3(X)  K.  Although  elastic  constants 
depend  on  temperature,  their  values  are  not  provided  over  a  wide  range  of  temperature.  Thus,  it  is 
likely  that  if  temperature-dependent  elastic  constants  become  available,  die  accuracy  of  calculation  can 
be  enhanced  particularly  at  low  temperatures.  At  the  same  time,  inclusitm  of  other  decay  mechanisms 
(i.e.,  other  dian  an  LO  phonon  decaying  into  two  LA  phonons)  reduces  the  lifetime  by  approximatdy 
0.5  ps  and,  thus,  furdier  improves  the  agreement  From  the  figure,  it  is  also  found  that  the  jdionon 
lifetime  does  not  depend  strongly  on  the  [dionon  wavevector.  As  a  result  the  k=0  value  provides  a 
good  approximation  for  the  LO  {dionon  lifetime  in  bulk  materials. 

Rionon-assisted  F-X  transition  rates  in  tvpe-II  supeiiattices 

The  relaxation  rates  of  (diotoexcited  electrons  is  an  important  issue  because  of  its  fiindamaital 
nature  as  wdl  as  for  possiUe  device  triplications.  Contrary  to  the  case  of  type-I  heterostructures,  die 
electrons  in  type-n  structures  are  excited  to  a  direct-gap  energy  level  (T)  in  one  layer  aid  then  can 
relax  to  an  indirea-gap  energy  level  (X)  in  the  adjacent  layer.  This  case  is  potraitially  significant  for 
device  apidications  since  electron  relaxation  occurs  in  conjunction  widi  real-space  transfer.  This  pro¬ 
cess,  which  is  normally  fotbidden,  is  now  possible  due  to  the  mixing  of  T  and  X  states  by  die  superlat- 
tice  potential  and  relaxation  of  momentum  conservation  due  to  interface  disorder  (lateral  mixing). 


Several  experimental  results  related  to  F-X  transfer  are  available  in  the  literature.  In  spite  of  the  spa¬ 
tial  charge  transfer,  the  relaxation  rate  can  be  very  high.  To  this  date,  there  has  been  relatively  little 
theoretical  efibrt  to  analyze  the  F-X  transfer  in  type-n  structures.  Recently,  we  have  calculated  the  F- 
X  scattering  rates  due  to  optical  phomn  emission  using  a  realistic  band  structure  model  (tight  binding), 
and  have  demonstrated  that  diis  is  an  important  mechanism  for  the  F-X  relaxation  process  in  type-II 
siqwilattices. 

The  electron-fAionon  interaction  in  low-dimensional  systems  is  altered  strongly  due  to  the 
confinement  of  the  carriers  and  the  confinement  of  the  phonms.  Usually,  a  single-band,  spherical 
effective-mass  model  is  used  for  the  description  of  the  confined  carrier  states.  It  is  commonly  assumed 
that  eadi  of  these  states  is  derived  from  bulk  states  of  a  givoi  symmetry  only  (e.g.,  F,  X,  or  L),  and 
levels  derived  from  different  bulk  states  do  not  interact  widi  each  other.  In  a  supeilattice,  momentum 
in  the  growth  direction  is  not  conserved  due  to  the  discontinuities  in  die  superiattice  potential  in  dns 
directkm.  Therefore,  bulk  states  of  the  same  energy  but  different  symmetry  can  craiide  with  each 
other.  As  a  result,  the  envelope-functicHi  approximation,  which  woiks  very  well  in  many  cases,  is  not 
suitable  in  this  case.  In  our  study,  a  second  nearest  neighbor  empirical  tight  binding  method  widi  an 
sp^  basis  was  used  for  describing  the  electronic  band  structure  [14].  The  mixing  of  the  F  and  X  val¬ 
leys  is  intrinsically  included  in  this  model  since  the  complete  band-structure  is  described  in  die  fight 
binding  n^hod.  For  a  proper  calculation  of  the  electron-phonon  scattering  rates,  the  effects  of  phonon 
confinement  should  also  be  taken  into  account  Although  we  have  developed  a  microscopic  model  for 
lattice  dynamics,  using  dus  model  iii  the  calculation  of  electron-fdionon  interactions  coupled  with  a 
resdisfic  etectronic  band  structure  requires  intensive  computafioa  Therefore,  for  die  description  of  [dio- 
non  OHifinement  we  used  the  dielectric  continuum  (slab)  model  because  of  its  simplicity  and  reason¬ 
able  accuracy.  In  the  didectric  ctmtinuum  model  picture,  the  confined  modes  are  totally  confined 
widiin  an  individual  layer,  and  die  fdionon  potential  of  one  type  of  layer  is  identically  zero  in  die  adja¬ 
cent  layer.  Hoioe,  (dionon  potentials  in  each  wdl  do  not  interact  with  eadi  other,  and  omsequently. 
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confined  i^kxkmi  modes  in  a  superiatdce  are  the  same  as  the  ones  for  a  quantum  well  in  this  picture. 
The  inteiface  phontm  modes  in  a  superiattice.  on  the  other  hand,  are  modified  due  to  the  periodicity  of 
the  stnicture  and  die  oveilapinng  of  potentials  [IS].  In  diis  case,  die  idionon  potential  extoids  over  the 
whole  superiattice.  The  resulting  Hamiltonian  due  to  interacting  interface  modes  is  very  complicated. 
To  simidify  the  calculations,  we  coisidered  only  the  case  where  the  phtmon  wavevector  along  the 
growdi  directitm  (i.e..  the  wavevector  in  the  superiattice  miniband)  is  zero.  Based  on  diis  interaction 
Hamiltcmian,  we  calculated  the  election  scattering  rate  due  to  optical  {Aionon  emission  following 
Fermi’s  Golden  Rule.  To  perform  the  summation,  we  made  use  of  the  fact  that  the  phonon  potentials 
vary  slowly  in  the  unit  cell  compared  with  die  atomic  orbitals  which  are  strongly  localized  at  atomic 
locations.  The  sum  over  all  possible  finite  states  within  the  2-D  Brillouin  zone  was  evaluated  numeri¬ 
cally  [16]. 

The  specific  structure  of  interest  in  our  initial  calculations  was  a  superiattice  with  M  mmolayers 
of  GaAs  and  N  monolayers  of  AlAs.  The  energy  levels  at  die  miniband  minimum  obtained  horn  die 
tight  binding  calculation  as  a  function  of  AlAs  layer  thickness  are  shown  in  Rg.  4.  The  GaAs  layer 
thickness  was  kqit  constant  at  8  monolayers.  The  levels  are  labeled  as  F  and  X  following  the  effective 
mass  notation,  although  each  level  is  actually  a  combination  of  these  two  bulk  states.  Here,  die  super- 
lattice  states  are  labeled  after  the  dominant  of  the  two  bulk  states  forming  them.  As  can  be  seen,  the 
effect  of  changing  the  AlAs  layer  thidmess  is  to  make  the  energy  difference,  hence,  the  interaction 
between  die  lowest  F  level  (Fi)  and  the  X  levels  vary.  Ihe  parity  behavior  of  states  needs  to  be  moi- 
ticmed  briefiy  because  of  the  effect  of  parity  on  the  overlap  of  electronic  wave  functions.  In  the  super- 
lattice,  die  electnmic  wave  functions  have  definite  parity  only  when  the  electrtm  is  at  the  miniband 
minimum  or  miniband  maximum.  When  the  initial  state  for  the  electron  is  assumed  to  be  the  bottom 
of  the  miniband,  the  parity  of  the  F|  state  is  even,  regardless  of  AlAs  layer  thickness.  The  parity  of 
the  X  levels,  on  the  other  hand,  dqpends  on  the  number  of  AlAs  monolayers.  The  X^  levels  with  odd 
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Hguie  4.  Energy  levels  obtained  from  the  tight  binding  calculatitm  as  a  fiinction  of  AlAs  layer  thick¬ 
ness.  The  GaAs  thickness  is  fixed  at  8  monolayers.  The  discrete  energy  tevels  are  maiVaH  as  F  and  X 
after  the  more  dominant  one  of  die  two  bulk  states  making  up  the  confined  stale. 


Figure  5.  Fj-Xi  tnmsition  rate  due  to  optical  phonon  emission  as  a  fiinctitm  of  AlAs  layer  thicki^s. 
The  GaAs  thideness  is  8  monolayers  throughout. 


(even)  i  have  the  same  (opposite)  parity  as  N.  The  parity  behavior  of  states  is  importam  because  the 
states  of  opposite  parity  do  not  mix  strongly. 

Scattering  due  to  the  onission  of  differait  phonon  modes  are  shown  in  Fig.  S.  The  results 
IKesented  are  the  sums  of  contributions  by  all  normal  modes  for  a  given  type  of  i^non.  In  calculat¬ 
ing  the  inteisubband  scattering  rates,  die  election  is  taken  to  be  initially  at  the  bottom  of  the  F)  sub¬ 
band.  It  then  transfers  to  the  X]  subband  by  emitting  an  optical  phonrai.  As  is  apparent  from  Figs.  4 
and  5,  die  scattering  rate  is  strongest  when  the  F  and  X  levels  are  closest  in  energy.  Under  this  crnidi- 
don,  the  levels  interact  very  strongly  with  each  other  and  the  F  level  has  large  X  character.  Iherefore, 
the  oveihq)  between  the  inidal  and  final  states  is  large.  This  interaction  and  the  resulting  large  overiap 
decreases  rapidly  as  die  energy  separation  increases.  This  is  one  of  die  two  reasons  for  strong  diidc- 
ness  dependence  (the  other  being  parity)  of  the  scaring  rate.  Among  different  phonon  modes,  die 
ALAs  confined  modes  are  strongest,  followed  by  die  AlAs  interface  modes.  Hie  GaAs  modes  are 
weaker,  witii  die  GaAs  confined  modes  being  the  weakest  This  is  due  to  the  fact  that  the  final  state 
X}  wave  function  is  strongly  confined  within  the  AlAs  layers.  Consequently,  the  overly  involving 
GaAs  confined  modes  is  very  small.  An  interesting  point  to  note  about  Fig.  S  is  that  the  dominant 
phonon  modes  alternate  between  being  even  and  odd.  This  is  because  die  parity  of  the  X  level  alter¬ 
nates  as  the  AlAs  thickness  is  changed  by  one  monolayer  each  time.  For  even  (odd)  AlAs  layer  diick- 
ness,  the  parity  of  the  lowest  X  level  is  even  (odd);  hoice,  the  even  (odd)  [dionon  modes  are  straig. 
The  effect  of  parity  can  be  also  be  observed  at  cross-over  points  where  the  F  and  X  level  oieigies 
bectnne  very  dose.  From  Fig.  4,  it  is  seen  that  die  two  levels  have  almost  die  same  energy  for  N=S, 
but  tiiey  do  not  mix  and  repel  each  other.  This  is  because  F]  and  X]  have  the  (^tposite  parities  for  this 
thidaiess.  For  other  diicknesses  where  the  two  levels  become  close,  however,  diey  do  iqiel  eadi  other 
because  they  have  die  same  parity.  This  is  the  reason  why  the  scattering  rate  for  Nsl9  is  mudi 
stronger  than  for  N=20.  Finally,  we  compared  our  results  with  an  eiqieriment  by  de  Paula  [17]  where 


20 


Ok  r  to  X  traosition  rate  via  i^ionon  emission  was  obtained  a  time-resolved  ami-Stokes  measure- 
menL  For  the  case  of  M=:8  and  Nsl4,  they  found  the  transfer  time  to  be  around  1  ps,  which  is  in 
agreement  with  the  rate  calculated  in  this  study. 

Extension  of  .Kronig-Pennv  model  for  T-X  mixing  in  type-Il  superlattices 

As  discussed  in  the  previous  section,  an  accurate  description  of  band  structure  including  T-X 
mixing  is  essential  for  the  investigation  of  electronic  and  optoelectronic  properties  in  type-II  superiat- 
tices.  For  the  conduction  bands  of  supeilattices  made  of  direct-gap  materials  (i.e.,  type-I),  the  Kronig- 
Poiney  has  been  used  satisfactorily  to  understand  many  important  proUems.  However,  other  compli¬ 
cated  m^hods  (such  as  the  empirical  tight  binding  method  used  in  the  previous  section)  are  required  to 
model  the  superiatdce  when  this  is  not  the  case.  An  accurate  yet  simple  model  for  treating  mixing 
effects  is  useful  to  make  theoretical  predictions  and  explain  the  experimental  results.  Therefore,  it  is 
very  desirable  to  extend  the  Kronig-Penney  model  so  that  mixings  of  different  valleys  may  be  treated. 
We  devdoped  such  an  extoision  in  whidi  the  energy  bands  and  wave  functions  in  type-II  superlattices 
are  calculated  adopting  a  simplified  tqiproach  formulated  by  Cho  and  Prucnal  [18]  altmg  with  the  boun¬ 
dary  amditions  introduced  by  Uu  [19]  in  treating  single  barrier  tunneling.  Our  model  has  die  merit  of 
being  easy  to  implement  and  permits  the  wave  functions  to  be  expressed  analytically.  Therefore,  it  can 
be  a  valuable  tool  for  simple  and  efficient  calculatkms  of  mixing  effects. 

In  a  superiatdce,  the  electronic  state  is  frequoidy  apinoximated  by  a  mixture  of  the  related  ztme- 
center  (T)  and  zcme-edge  (X)  bulk  states.  Accordingly,  the  wave  foncdon,  y,  is  written  as  die  sum  of 
the  products  of  slowly  varying  envelope  foncdons  and  band  edge  periodic  foncdons  (Bloch  funcdons): 
i.e.,  y  =  ypUr  +  yxux>  where  ur  and  ux  are  die  Bloch  funcdons  for  die  F  and  X  valleys,  respecdvely. 
In  vector  notation,  the  wave  function  may  be  exixessed  as  y  =  (yr.yx)-  The  coupling  of  F  and  X  val¬ 
leys  due  to  the  InOken  translational  symmetry  of  the  crystal  can  be  modeled  by  die  off-diagmial  ele¬ 
ments  in  the  potoitial  energy  [19].  For  simidicity,  these  off-diagonal  terms  are  refnesented  by  ddta 
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functkms  centered  at  the  interfaces: 


V  = 


VrCz)  aS(z) 

p£(z)  Vx(2). 


(13) 


where  z  is  the  directiai  of  crystal  growth,  a  is  a  parameter  quantifying  mixing,  S(z)  is  the  Dirac  delta 
function,  and  Vr(z)  and  Vx(z)  define  the  F-potential  discontinuity  and  X-potendal  discontinuity, 
respectively.  The  wave  (envelope)  function  continuity  condition  (at  z=0)  is. 
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Vx 

0- 
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(14) 


The  other  boundary  condition  can  be  derived  by  integration  of  the  Schrbdinger  equation  across  the 
interface: 


1  ayr  I 

2  mr(z)  dz  ^ 
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+  cniTx  =  0  , 
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+  ayr  =  0  . 


(15) 


As  can  be  seen,  introduction  of  the  delta  function  potential  results  in  a  boundary  condition  in  which  the 
r  and  X  vaUeys  are  coupled  across  the  interface.  In  obtaining  this  result,  it  is  assumed  that  the  Bloch 
fimctkms  of  the  respective  band  edges  in  the  two  constituent  materials  do  not  differ  significantly. 


With  our  model,  we  obtained  arudytical  expressions  for  the  envelop  functions  and  a  transcenden¬ 
tal  equatkm  for  die  energy  levels.  The  energies  of  die  lowest  F  and  X  levels  are  shown  in  Fig.  6  as  a 
function  of  layer  diickness  with  and  without  mixing.  In  this  figure,  a  equals  0.1  eV-A;  as  discussed 
below,  this  is  a  reasonable  value  in  light  of  the  favorable  agreement  with  numerous  experimental  and 
dieotedcal  results.  As  can  be  noticed,  the  effect  of  mixing  is  to  create  an  anticrossing  behavior  where 
dw  two  levels  are  expected  to  meet.  This  important  feature  is  cortectiy  reflected  in  our  treatment  We 
also  demonstrated  that  die  wave  functions  in  type-n  superlattices  obtained  fnmi  conventional  Kroriig- 
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Figure  6.  Bneigies  at  the  miniband  minimum  for  the  two  lowest  levels  in  a  GaAs/AlAs  superiattice 
as  a  fiinctKMi  of  layer  Oiickness.  Hie  GaAs  and  AlAs  layers  are  assumed  to  be  of  equal  diickness. 
Ifigher  energy  levels  are  omitted  in  dds  figure  for  darity.  Solid  lines  show  the  bands  with  a=0.1  eV- 
A,  udiile  die  energies  for  dashed  lines  are  obtained  by  setting  a  equal  to  0  (i.e..  conventitmal  Kronig- 
Pmney  model).  The  inset  shows  the  anticrossing  behavior  in  detail  and  has  the  same  labds  and  the 
units  ra  die  main  figure. 
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Penney  modd  are  considerably  in  enor  when  the  layer  thicknesses  are  such  that  strong  mixing  of  F 
and  X  states  occurs.  Near  the  point  of  andcrossing.  it  is  imperative  that  the  interaction  of  the  two  val¬ 
leys  be  considered  to  have  correct  wave  functions.  Our  model,  as  opposed  to  the  conventional 
Kronig-Penney  model,  yields  af^ropriate  wave  functions  for  the  mixed  states.  We  obtained  the  value 
of  our  mixing  parameter  (the  only  adjustable  parameter  in  our  model)  by  comparison  with  bodi  experi¬ 
ments  and  more  soi^sticated  theoretical  models.  The  values  of  a  deduced  from  a  variety  of  sources 
conveiged  around  a  single  value  (0.1  eV-A)  of  diis  parameter.  As  a  result,  we  showed  that  the  simile 
and  efficient  model  we  developed  could  reasonably  reproduce  the  impoitant  aspects  of  valley  mixing 
effects  in  Q^-II  superiattices. 

Theory  of  Bloch  electron  dynamics  in  electric  fields 

Bloch  electron  dynamics  in  electric  fileds  has  been  a  subject  of  great  interest  dating  back  to  the 
eailiest  ai^cations  of  quantum  mechanics  to  solid  state  i^ysics  [1].  Even  more  recenfly,  as  modem 
fabricatkm  technologies  continue  to  drive  the  study  of  solid  state  transport  into  the  nanometer  domain, 
many  new  and  interesting  questitHis  have  emerged  concerning  die  solid  state  dynamics  and  quantum 
transport  of  carriers  in  band-engineered  superiattices  and  tailored  periodic  solids. 

In  die  past  year,  we  have  developed  a  novel  formalism  for  treating  Bloch  electron  dynamics  and 
quantum  transport  in  a  superimposed  uniform  and  oscillatory  electric  field  ]^t)  =  ^  +  ^iCOScoL  In  die 
formalism,  the  electric  field  is  described  throu^  the  use  of  the  vector  potential,  and  a  basis  set  of 
localized,  electric  field-dependent  Wannier,  and  related  envelope  function  is  developed  and  utilized  to 
describe  the  dynamics.  From  the  dieory,  a  complete  analysis  of  single-band  and  muld-band  processes 
is  derived  for  a  Blodi  electron  in  the  presence  of  a  superimposed  uniform  and  oscillatory  electric  field, 
hi  the  single-band  approximation,  it  is  shown  that  the  tuning  of  the  oscillatory  field  fiequeircy  <o  to 
some  multiple  of  the  Blodi  frequency  cog  »  eEoa/h  results  in  die  resonant  displacement  of  the  elec¬ 
tron  fipom  a  ^en  site  to  die  tuned  muldple  of  latdce  dtes  along  with  a  concomitant 
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(tevelopment  of  a  DC  vdodly  in  the  direction  of  the  electric  field  for  the  initial  Bloch  state, 

^  Z  c„(inNoaXi)  (®M)  (-ir  Jm(®Noa)  sin(mNoK^a) .  (16) 

“  noO 

where  Jn  is  the  Bessel  function  of  integer  order  m.  o»E)a/tiO),  and  die  energy  hand  function  and  its 
Fourier  components  are  defined  as 

e.(it)  =  e.(K,JtD  =  £  c.(l.Jti)  e**^  =  £  e*^  y  c.'^  e^^^-L . 

1.  a 

On  die  odier  hand,  for  a  fixed  oscillatory  field  frequency  co,  the  DC  current  can  occur  in  the  direction 
of  the  electric  field  by  tuning  the  Bloch  frequency  0%  (DC  field  strength  E^)  to  some  multiple  of  the 
oscillaloty  field  frequency  (i.e.,  (a^MoCa).  The  DC  velocity  for  initial  Bloch  state  is  found  to  be 

Z  c.(l,.iU  Jy^Oxa)  lx  sina,K„a)  .  (18) 

“  IpO 

For  the  one-dimensional  nearest-nei^ibor  fight-binding  model,  the  DC  velocity  is 

^  j^(„)  sij,(K^)  ^  (19) 

whidi  resembles  the  usual  Bloch  velocity  with  a  deformed  band  parameter 

Further,  firom  the  multi-band  arudysis,  a  ^neralized  transition  rate  is  derived  in  terms  of  the  DC 
field,  oscillatory  field,  and  dK  relevant  band  parameters,  showing  all  possible  transitions  include  Zener 
tunneling,  photon  assisted  tunneling,  mulfiphoton  onission/absorption,  as  well  as  transitions  between 
stiUes  of  the  Wanmo’  laddo*.  Based  on  dre  methodology  whidi  was  devdoped  and  presented  in  our 
earlier  rqiott,  the  muM-band  coupling  is  treated  in  the  Wigner-Weisskopf  s^ipioximation,  whidi  allows 
for  the  analysis  of  die  Itmg-time,  time-dependent  tunneling  diaracteristics  of  an  electron  transition  out 
of  an  initially  occuffied  band  due  to  power  absoibed  by  the  electric  fidd,  while  preserving  conservation 
of  total  transition  probalfility  over  the  comfdete  set  of  exdted  bands.  The  total  transition  rate  for  an 
electron  to  tunnel  out  of  die  inifially  occupied  band  is  derived,  and  the  results  show  the  dqiendence  of 
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the  transitkm  rate  on  the  initial  ccHidition  of  the  system  and  the  decay  rate  of  the  initial  band  amplitude, 
YaOt),  whidi  is  given  as. 


.  n'#n  " 


4  n 


(20) 


i.^m  *1  Ua 


It  is  noted  that  Yn^)  consists  of  several  S  functions,  which  indicate  the  possible  physical  processes  for 
a  Bloch  etectnm  in  a  homogoieous,  time-vaiying  electric  field.  The  first  term  in  Ya(i^)  indicates  the 
inteiband  transition  due  to  Zener  tunneling  betweot  the  deformed  bartds  n  and  n^  which  is  proportional 

to  Here.  D„<l,J^i)=  lae"^'**'.  €„(qjti)=  e„<q4tx).  and  6^'=  Mit)  -  ^(it)l/lr  with  the 

4 


defonned  enetgy  band  function. 


l,m 


(21) 


where  IxODb  mm  =  0.  The  second  and  third  terms  in  Yn(^)  indicate  the  tranddons  between  the  two 
deformed  bands  dirough  absoiptioiVemission  of  a  photon,  wdiich  are  proportitmal  to  Ef,  and  foe  last 
term  suggests  the  inteiband  transitions  between  foe  deformed  bands  via  Watuiier-Staik  transitions,  as 
wen  as  multiphonon  absoiption/emission.  The  Franz-Keldyfo  effect  can  dso  be  explained  from  foe 
results. 


In  addition  to  Bloch  electron  dynamics  in  a  homogoieous,  time  varying  decliic  field,  we  have 
investigated  the  effea  of  localized  impurities  on  the  transport  of  an  election  in  a  uniform  electric  field 
Sq.  We  have  developed  a  theory  to  describe  Blodi  elecmxi  transport  forough  spatiaHy  localized  inho¬ 
mogeneous  potentid  barriers.  In  ^qfiying  foe  theory,  field-dependent  envelope  functions  and  currents 
for  a  Blodr  dlectrrm  tunnding  through  single  and  double  Slater-Koster  potential  barriers  are  derived 
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using  the  Wannier  representaticni  within  the  one-dimensional  nearest-neighbor  tight-binding  approxima¬ 
tion.  The  exaa  solutions  for  the  envelope  functions  in  the  Wannier  representation  are  expressed  in  an 
integral  form  through  the  use  of  the  Lr^lace  transform  solutions  for  the  envelope  functions  at  the 
impurity  sites.  Eaiiy  time  solutions  are  obtained  for  die  envelope  functions  and  the  velocity  is  calcu¬ 
lated  for  the  given  fields.  The  asymptotical  behaviors  for  hi^-field  and  low-field  limits  ate  observed. 

To  study  the  effects  of  the  localized  potoidal  barriers  on  die  Wannier-Stark  levels,  we  have 
derived  the  electric  field  dependent  Green's  function  for  the  one-dimensional  nearest-neighbor  tight- 
binding  band,  e(it)  =  +  2eiCOsKa,  which  is 


^ 


smicy 


(22) 


where  y  =  (&-eoVli<DB-  The  Stark  energy  spectra  for  single  and  double  Slater-Koster  impurities  are 
calculated  using  the  field-dependent  Green’s  fimction.  The  results  show  that  for  relatively  low  poten¬ 
tial  barriers  and  high  electric  field,  the  energy  levels  ate  almost  uniform,  resembling  the  Warmier-Stark 
ladder,  and  the  barriers  act  as  perturbations;  for  relatively  low  field  and  high  potential  barriers,  die 
Wannier-Stark  uniformity  of  the  levels  becomes  distorted  (see  Fig.  7).  Although  the  illustrations  cited 
here  ate  for  Slater-Koster  barriers,  the  calculation  is  easily  extended  to  general,  localized  potentials. 


2.3  Publications  and  Presentations 

During  the  last  year,  this  program  has  resulted  in  nine  refereed  publications  in  the  literature,  five 
additional  manuscripts  are  curtendy  in  press,  and  one  more  has  been  submitted  to  a  major  technical 
JoumaL  Ten  presoitations  and  invited  talks  have  been  given  at  conferences  and  workshops.  The  fol¬ 
lowing  paragraphs  summarize  the  publications  and  presentations  made  under  this  program  during  the 
past  contract  period. 
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Rguie  7.  Wannio’-Staik  eneigy  levels  for  two  impurities  of  two  lattice  ^Kicing  apart  with  die  teladve 
impurity  strength  of  (a)  V,/2ei  s  OJ  and  (b)  Vg/Tei  s  l.O. 
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An  eztenaioa  of  the  Kronig-Penney  model  to  treat  mixings  of  different  valleys  in  supeiiattices 
it  pfcaented.  The  value  (ff  the  mixing  parameter  a  used  in  the  model  is  deduct^  by  comparison 
with  moR  ao|dusticated  theory  and  various  experiments  for  the  OaAs/AlAs  mate^  system.  It 
is  shown  that  with  the  me  <ff  a  single  value  for  the  mixing  parameter,  this  simple  and  efficient 
model  can  reasonaUy  rqjtoduoe  the  inqMrtant  aspects  of  valley  mixing  effects  in  OaAs/AlAs 
siqMflattioes  over  a  wide  range  of  layer  thicknesses. 


The  GaAa/Al/3ai_;(As  superlattices  with  x>0.43 
have  been  studied  extensively  becinise  of  their  new  physical 
phenomena  and  possilde  apidkations.*  In  these  superlat¬ 
tices,  for  a  certain  range  of  layer  thicknesses,  F  and  X 
vaDqrs  can  coiqik  with  each  other  strongly  due  to  the 
broka  tranalatioiial  symmetry  of  the  host  lattice  caused  by 
the  potential  disoontmoity  at  the  interfaces.  Recently,  there 
has  been  coosideraMe  interest  in  this  subfect^'^  For  the 
conduction  bands  of  superlattioes  moAo  of  direct-g^>  ma¬ 
terials,  the  Kranig-Peimey  model  can  be  used  satutiGurto- 
rily.  However,  other  complicated  methods*"**  are  required 
to  modd  the  snperlattice  when  this  is  not  the  case.  An 
accurate  yet  simple  model  for  treating  mixing  effects  is 
uaeftal  to  make  theoretical  predictions  and  explain  the  ex¬ 
perimental  results.  Therefore,  it  is  very  desira^  to  extend 
the  Kronig-Pennqr  model  so  that  mixings  of  different  val- 
Iqrs  may  be  treated.  In  this  work,  we  develop  such  an 
extension,  in  which  the  energy  bands  and  mixed-symmetry 
wave  functions  in  superiattices  are  calculated  adopting  a 
simplifled  apfnoach  formulated  by  Cho  and  Prucnal,*^ 
alo^  with  the  boundary  conditkns  introduced  by  Liu^  in 
treating  single-barrier  tunneling.  It  is  our  aim  to  demon¬ 
strate  that  this  simple  model  with  a  single  value  (ff  the 
mixing  parameter  for  modding  the  strength  of  F-X  inter¬ 
action  can  reproduce  the  important  aqiects  of  the  experi- 
mentd  results  as  wdl  as  th^  at  the  mote  complicated 
theory  over  a  wide  range  of  dimenskmd  scales.  This  aim- 
lifted  modd  uses  the  effective-mass  theory,  neglects  the 
nonparabolicity  effects,  and  does  not  ded  with  lateral  mix¬ 
ing  due  to  intoffice  roughness**  however,  it  has  the  merit 
of  being  easy  to  imidement  and  permits  the  wave  functions 
to  be  exptei^  analytkally.  Therefore,  it  can  be  a  vduable 
tool  for  simfde  and  efficient  calculation  of  mixing  effects. 

In  a  superiattioe,  the  electron  state  is  a  mixture  of  the 
zone-center  F  and  zone-edge  related  bulk  states.  Accord¬ 
ingly,  the  wave  Amction  ^  is  written  as  the  sum  of  the 
IMoffiicts  of  slowly  varying  envelope  flinetions  and  band- 
edge  periodic  functions  (Bloch  fimctions):  i.e., 
where  Up  and  ux  are  the  Bloch  functions 


for  the  F  and  X  valleys,  respectivdy.  In  vector  notation, 
the  wave  fonctum  may  be  expressed  as  1°  * 

superiattice,  the  six-fold  degeneracy  (ff  ^  va^  is  re¬ 
moved,  with  the  formation  of  an  jf,  doublet  (momentum 
ahmg  the  growth  axis),  and  an  Xg^  quadruplet  (momen¬ 
tum  in  the  plane).  The  Xg^  minima  are  not  coupled  to  the 
F  valley  due  to  the  conaervatkm  of  lateral  momentum. 
However,  the  F  and  Xg  valleys  ate  coupled  with  each  other 
since  momentum  in  the  z  direction  is  not  conserved  be¬ 
cause  of  the  discontinuities  in  the  superiattice  potentiaL 
Accordingly,  out  of  all  the  AT  states,  o^  the  Xg  states  ate 
mixed  with  the  F  states  by  the  perfect  superiattice  poten¬ 
tial,  and,  herein,  the  X  point  will  be  used  to  mean  Xg  par¬ 
ticularly.  The  coupling  of  F  and  X  valleys  due  to  the  bro¬ 
ken  translatioDal  symmetry  of  the  crystal  can  be  modded 
by  the  off-diagonal  dements  in  the  potential  energy.^  For 
simplidty,  these  terms  are  reinesented  by  ddta  functions 
centered  at  the  interfaces. 


/Kr(z)  a6(z)\ 
^o6(z)  Kx(z)j' 


(1) 


where  a  is  a  parameter  quantifying  mixing,  ^(x)  is  the 
Dirac  delta  fimctimi,  and  Vpiz)  and  Kx(z)  de&ie  the 
F -potential  discontinuity  and  JF-potentid  discontinuity,  re¬ 
spectivdy.  The  wave  (envekqwj-function  continuity  om- 
dition  (at  z=0)  is  (tprf  'l>x)o+-  The  otJ»er 

boundary  condition  can  be  derived  by  integration  of  the 
Schrodinger  equation  across  the  interface. 


1^1  1 
2  ^mr(z)  dx 


0+ 


^  /  1  ^zj 

2  ^mx(s)  ar 


1  d^i  \ 
1  dtjixx  \ 


(2) 


As  can  be  seen,  introductum  of  the  ddta-fonction  potential 
results  in  a  boundary  comlititm  in  which  the  F  and  X 
valleys  are  coujded  across  the  mterfoce.  In  obtaining  this 
result,  it  is  assumed  that  the  Bloch  functions  of  the  respec- 
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A  micfoecoiric  model  of  optical  phonons  in  dimensionally  confined  structures  is  formulated 

and  an^ed  to  calculate  electnM-optiimd'phonon  scatterini  rates  in  OaAs/AlAs  quantum  wells.  For 
this  simpltaert  model  which  circumvents  performing  a  complicated  ab  initio  calculation  of  the  force  con¬ 
stants  at  the  interfiice,  it  is  demonstrated  that  the  resulting  dispersion  relation  and  scattoing  rates  for 
electron-optkal-phonon  interactions  agree  very  well  with  those  obtained  from  detailed  ab  initio  studies. 

It  is  also  shown  that  for  OaAs/AlAs  structures,  the  macroscopic  dielectric  continuum  model  provides  a 
good  ^tprozimation  to  the  scattering  rate  predicted  by  the  microscopic  models. 


The  electroiuc  and  optical  properties  of  semiconductor 
saperlattkes  (SL*s)  and  quantum  wells  (QW’s)  have  been 
investigated  extensively.  A  principal  advantage  of  using 
such  heterostructures  results  from  the  ability  to  tailor  the 
electronic  and  optical  properties  of  the  structures  for 
realizing  a  potentially  vast  array  of  high-performance 
eleetrmic  a^  optodectronie  devices.  To  ftilly  under* 
stand  and  utilize  the  properties  of  these  nanometer'4cale 
heterostructures,  it  is  necessary  to  develop  formalisms  for 
studying  confinement  effects  as  well  as  picosecond  and 
sobpkoeecond  processes.  It  has  been  known  for  many 
years  that  the  scattering  by  polar-optical-phonon  modes 
is  an  impwtant  energy-loss  mechanto  for  electrons  in  a 
wide  variety  of  III-V  semiconductor  devices.  However, 
effects  of  confinement  on  these  phonon  modes  have  been 
investigated  extensively  only  in  Ae  past  several  years. 

In  recent  years,  a  numb^  of  models  has  been  put  for¬ 
ward  to  explain  electron-optical-phonon  interactions  in 
reduced  dimensional  systems.  Tliey  can  be  broadly 
classified  in  two  categories:  macroscopic' and  micro¬ 
scopic. Macroscopic  models  ignore  the  effect  of  in¬ 
dividual  layers  of  atoms  but  they  have  the  considerable 
advantage  of  making  the  interaction  calculation  very  sim¬ 
ple.  Among  these  macroscopic  models  are  the  dielectric 
continuum  model*''’  (slab  nKxlel,  which  uses  purely  elec¬ 
trostatic  boundary  conditions),  hydrodynamic  model,^ 
hybrid  model,*  and  a  recent  disperwve  continuum  treat¬ 
ment  of  Nash.’  In  some  parameter  regimes,  these  models 
are  Burly  accurate  and  provide  good  estimates  of  energy- 
loss  rates.  However,  scaling  of  the  electron-optical- 
phonon  interaction  with  diminiihing  device  length 
presentt  a  serious  challenge  to  the  accurate  use  of  such 
models.  As  a  result,  there  has  recently  been  an  increasing 
need  for  more  rigorous  analysb  and  detailed  knowledge 
of  electron-optical-phonon  interactions  in  reduced  di¬ 
mensional  systems.  This  has  been  the  main  motivation 
for  the  emergence  of  ab  initio  microscopic  models.'^ 
Thoui^  such  nmdels  provide  the  most  accurate  analysis 


of  the  structure,  they  have  not  been  used  extensively. 
This  can  be  attributed  to  the  fimt  that  the  ab  initio  micro¬ 
scopic  analysis  involves  very  arduous  aitd  time  consum¬ 
ing  first-principle  calculations  of  lattice  d)mamics'*''* 
rather  than  employing  adjustable  parameters.'’~  *’ 

Precise  ab  initio  calculations  of  force  constants  at  the 
interface  may  not  be  essential  for  moot  (rfthe  heterostruc- 
tures  except  those  involving  extremely  thin  layers.  It  is 
well  known  that  even  a  simple  linear-chain  model  with 
nearest-neighbor  force  constants  can  predict  the  zone- 
center  LO-phonon  frequencies  in  a  SL  with  a  reasonable 
accuracy  except  in  the  cases  where  layers  are  single 
monolayer  thick.  Such  an  approximate  model  is  based  on 
the  assumption  that  atomic-force  constants  at  heterojunc¬ 
tion  interfaces  are  identical  to  those  of  the  bulk  or  of  uni¬ 
form  pseudomorphic  layers.  In  a  qualitative  analysis  of 
the  effect  of  varying  force  constants  at  the  hetmjunction 
interfaces  of  a  strained  lajrer,  short-period,  GaAs/GaP 
SL  with  two  monolayers  per  SL  layer,  it  was  also  found 
that  frequencies  of  the  confined  phonon  modes  are  cmly 
weakly  dependent  on  the  variations  in  the  interfaciid 
force  constants."  The  variation  in  interfacial  force  con¬ 
stants  by  values  as  extreme  as  10%  results  in  less  than 
about  8  2%  change  in  the  frequencies  of  confined  phonon 
modes.  It  should  be  noted  that  u  a  practical  matter, 
changes  in  the  frequencies  of  the  confined  LO-phonon 
modes  will  be  considerably  less  than  2%,  dnee  in  most 
SL’s  and  QW’s  the  ratio  of  the  number  of  bonds  at  the  in¬ 
terfaces  to  the  number  of  bonds  one  or  more  monolayers 
away  from  the  interfaces  is  less  than  that  for  the  case 
where  each  layer  is  two  monolayers  thick. 

Based  on  the  results  of  Ref.  18,  as  well  as  on  support¬ 
ing  observations  from  other  invest^tors,"  we  have  for¬ 
mulated  a  simplified  microscopic  model  which  facilitates 
the  accurate  modeling  of  confoied  and  interface  phonons 
without  ab  initio  calculations  of  force  constanu.  The 
valence-shell  model  developed  by  Kune  and  Nielson  for 
bulk^  has  been  extended  for  the  SL/QW  structures.  In- 
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Many  proposed  applications  of  mesoscqiic  elecmmic  structures  involve  carrier 
transport  at  low  temperatures  and  low  carrier  energies;  fiequently,  the  regime  of  interest 
IS  one  where  dimendonal  confinement  modifies  the  phase  spMct  substantially.  In  this  low 
temperature,  low  energy  reginie  [1-6].  acoustic  phomms  play  an  enhanced  in  carrier 
scattering  ioA  may  dominate  over  the  scattering  of  carriers  by  optical  phonons, 
pyrthennoie.  in  nanoscale  structures  it  is  possible  that  phase  space  restrictioos  may 
weaken  or  fortiid  optical  phonon  scattering  processes  tiua  woiud  normally  dominate  in 
bulk  structures.  In  recent  j«an.  there  has  been  an  extensive  literature  on  the  role  of 
dimensional  confinement  m  modifying  optical  phonon  modes  and  their  interactions  with 
durge  canien  (see,  for  example.  Refs.  7-1 1  arid  the  numerous  pliers  references  therein). 
Hoe^ver,  there  are  relatively  few  treatments  dealing  with  the  role  of  dimensional 
confinement  in  modifying  acoustic  phonon  modes  and  their  interactions  with  charge 
catTkn  [2-4].  While  dime  is  an  extensive  literature  tm  the  theory  of  acoustic  modn  in 
conventional  waveguides,  resonators  and  related  structures  [5].  no  efforts  have  been 
reported  to  formulate  a  theory  (rf  acoustic  phonons  in  nanos^  structures  where  both 
roontm  confinement  and  a  quantum  mechanical  tteannent  of  phonon  normalization  ate 
both  essential;  Gonstantinou  has,  however,  discussed  the  unnonnalized  acoustic  phonon 
modes  in  cylindrical  polar  aemicondoctor  quannim  wires  [12]. 

In  recent  ytm,  considerable  theoretical  and  experimental  effort  has  focused  oa 
understanding  the  role  of  confined  and  interface  longitudinal-optical  (LO)  phtmons  in 
nanoscale  and  mesoscopic  semiconductor  strucnires  [13].  Efowever,  relatively  little 
attention  has  been  given  to  theoretical  treatments  of  acoustic  phonon  confinement  or  to 
die  interaction  of  confined  acoustic  phonons  with  charge  catnen.  The  need  for  sudi 
theoretical  treatments  has  been  underscored  recently  by  experimental  studies  providing 
both  direct  and  indirect  [3,4]  evidence  of  the  importance  or  acoustic  phonon  confinement 
in  reduced  dimmisional  elMtronic  structures. 

In  this  paper,  by  appropriately  quantizing  the  acoustic  phonon  diqilacements  we 
have  obtained  the  contectly  nonnalized  expressions  for  acoustic  phonon  confined  in  a 
number  of  dimensionally  confined  nanostructures.  As  for  the  case  trfLO  phonon  modes 
in  rectangular  quantum  wires  and  quantum  dots,  there  are  no  exact  sedutions  for  die 
complete  set  <»  phmioo  modes;  however,  as  for  the  case  of  LO  phonon  modes,  the 
qrprmdmate  modes  presented  in  this  work  provide  simple  and  useful  mquessions  which 
are  well  suited  for  modeling  the  interaction  of  carriers  with  acoustic  phonons.  To 
investigate  the  e^o  of  reduced  dimensionality  on  the  coupling  between  acoustic 
phormns  and  chu^  carrien,  we  have  formulate  the  interaction  Hamiltonians  for  bodi 
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NOTE 

GENERALIZED  PIEZOELECTRIC  SCATTERING  RATE  FOR 
ELECTRONS  IN  A  TWO-DIME?>»10NAL  ELECTRON  GAS 

(XaMtocri  2t  May  1993;  M  rtetiedform  l9Jme  1993) 


mnooucTiON 

Many  envMoiMd  applKaboos  of  naMaeofMC  MoaiooiidiKtor 
muctura  aic  bawd  on  the  iniMpon  of  charpe  carrien  at 
low  tamperatufM  and  at  low  carrier  enerfies.  At  mch  low 
leaaperatuKi  and  b>w  carrier  enerpee,  the  aeatiering  of 
carrien  firom  acouetic  phonona  nay  doninalc  over  carrier 
eeatlcfiiig  from  optical  plio«oiM(l-^  Traniport  of  carrien 
in  a  two-dhaenrional  eleciran  gat  ia  of  ipe^l  importance 
within  the  field  of  carrier  tranaport  in  aemiconductora.  In 
thia  connection,  Prioe(I]  haa  given  an  cacettent  comprehen- 
aive  ticabnent  of  the  of  phonon  acattering  ia  twosli- 
I  dectroa  tranaport  in  polar  aemioondiictor  layers; 
apecMically.  for  the  case  of  a  twoKlimenaional  electron  gas, 
Mae  has  given  a  detailed  account  of  etectrott  acattering  by 
the  Frahh^  deformation  and  piamelectric  potentials.  In 
the  case  ofpmirwlaclric  acattering  of  electrons  from  aoouatic 
phonona.  Price  has  used  the  pieioeleetric  Constanta  givca  by 
Zook(2)  in  the  form  aummatiaed  by  Heann^3].  As  dit- 
caaaed  by  Price,  the  pieaodectric  oonstantt  gricn  by  Zodc 
were  derived  under  the  assumption  that  tlw  daatic  an- 
iaocropy  of  the  crystal  it  amaO.  Price  has  made  the  adriitional 
approitiaaation  of  awaragiag  Zook’s  pieaoelectric  conatanta 
over  the  aamuthal  dhectiont  in  the  fdane  of  the  two-dhnen- 
aiooal  dectroa  gaa  In  tUa  work.  Price’s  treatment  of 
pieaoeiecttie  acattering  is  geaerafiaed  by  not  taking  such  an 
average  over  thsae  azimiMhal  directkma. 


DBOUaSION 

As  diacumed  previoualy(l],  the  pieaoeiecttie  acattering  rate 
function,  i9'(l,2),  for  dectrona  in  a  twosiimeaaioiial 
electron  gu  interacting  with  the  acoustic  phonon  in  a 
aiacblende  crystal  may  be  written  u: 


»'(l.2)-y^2d(£(l)-£(2)J. 

(I) 

where: 

with: 

2  gVp»?  P*f/ 

(2) 

(3«) 

(3b) 

and 

^-36 

'  (q’  +  gV 

(4«) 

2^. + >«.  “  + eje? + Q-9^ 

(4b) 

in  these  reauitt  IPfl,})  equals  the  rale  of  trandtions  from 
initial  state.  I,  to  final  state,  2,  per  unit  volume  of  k  apace, 
g  is  the  acoustic  phonon  wavcvecior  normal  to  a  two-dhnen- 
aionat  dactron  gas  in  the  (100)  plane.  0,  aaO  (2-  are  the 
phonon  wavevacton  in  the  (100)  plane.  k  is 


Boltanann’a  constant,  r  is  the  temperature  of  the  lattice,  e 
is  the  charge  on  an  electron.  A, 4  is  tte  piezoelectric  constant 
which  has  a  value  of  roughly  i.S7  x  10*  V/m  for  CaAs,  p  is 
the  mass  density,  and  i,(r,)  is  the  longitudinal  (transverae) 
velocity  of  sound  in  tte  lattice.  The  limitaliotts  of 
eqos  (l>-(4b)  have  been  discussed  previously  in  Ref.  [I] 
where  it  is  noted  that  a  linear  proportionality  belwen 
phonon  frequency  and  phonon  wavevector  has  been  as¬ 
sumed.  Upon  performing  the  indicated  integrations  over  q: 


and 


**  g?g 

4  (Q)  +  Q\^' 


(5a) 


la. OrQ- 
e  ^  2  (Qf  e?)w  ■  4  (c;  ' 

Thus: 

kr(ek„)>f/9_  QIQ]  \\ 

i  L\* 

/a  I  3»  g;g»  \|-| 

V2(C?+e?)'" 


As  in  Refs  [1-3],  the  effects  of  screening  have  not  been 
considered  in  deriving  eqn  (6). 

In  the  limit  where  the  average  is  taken  over  azimuthal 


directions(l]: 

g?g?-g*/8. 

(7a) 

and 

g?-g?-gV2. 

(7b) 

with  the  result  that: 

kneh„fn(9  1  13  1  \ 

2  gV32pr?  32p*f/ 

(8) 

as  miginally  derived  by  Price. 

In  deriving  the  generalized  piezoelectric  scattering  rate  for 
a  two-dimensional  electron  gas,  it  has  been  assumed  that  the 
stiffness  constants  for  the  zincblende  crystal  are  isotropic. 
Indeed,  for  each  of  CaAs,  GaP,  InSb,  InAs  and  InP  the 
three  independent  experimentally-determined  stiffness  ctm- 
stants  difl^  in  magnitude  by  le»  than  3S%  from  stiffness 
constants  consistent  with  an  isotropic  medium.  On  the  other 
ha^  the  generalimd  expression  for  may  vary  from 
zero  to  twice  its  azimutlwlly  averaged  value,  9itl32Q  de¬ 
pending  upon  the  relative  magnitudes  of  g,  and  g..  While 
applications  such  as  calculations  of  m^lity  not 
considered  in  this  work,  the  fact  that  the  generalized  ex¬ 
pression  for  xA/g  may  vanish  as  g,  or  g.  goes  to  zero  may 
be  of  special  significance  of  the  frequentiy-discuased  cases 
where  tlK  dominant  charge  carrier  momentum  in  the  two-di¬ 
mensional  electron  gas  or  the  phonon  momentum  is  highly 
anisotropic;  in  particular,  mobilitim  for  such  systems  may 
be  enhanced  since  the  longitudinal  scattering  terms  may  be 
suppressed  as  a  result  of  the  anisotropic  form  of  eqn  (5a). 
TM  variation  in  tBJQ  is  mudi  less  pronounced  having  a 
range  from  I0n/32g  to  lfiii/32g  dep^ing  upon  the  rela¬ 
tive  magnitudes  of  g,  and  g.. 
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Multiband  theory  of  Bloch  electron 
dynamics  In  electric  fields 
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tU8  Army  Raaoareh  Offloo,  Raaearch  Triangle  Park,  NC  2770S-2211,  USA 


Afealraei  A  novel  multiband  theory  of  Bloch  electron  dynamica  in  homogeneous 
electric  flelda  of  arbitrary  atrength  and  time  dependence  is  presented.  In  this 
formaltam,  the  electric  field  is  described  through  the  use  of  the  vector  potential. 
Multibend  coupling  ia  treated  through  the  uae  of  the  Wigner-Weiaakopf 
approximation,  thus  allowing  for  a  Bloch  electron  transition  out  of  the  initial  band 
due  to  the  power  afaeorbed  by  the  electric  field:  also,  the  approximation  ensures 
conservation  of  the  total  transition  probability  over  the  complete  set  of  excited 
bands.  The  choice  of  the  vector  potential  gauge  leads  to  a  natural  set  of  extended 
time  dependent  basis  functions  tor  describing  Bloch  electron  dynamics  in  a 
homogeneous  electric  field;  an  associated  basis  set  of  localized, 
electrio-field-dependent  Wannier  and  related  envelope  functions  are  developed 
and  utilized  in  the  analysis  to  demonstrate  the  inherent  localization  manifest  in 
Btoch  dynamlea  in  the  presence  of  relatively  strong  electric  fields.  From  the 
ttwory,  a  generalized  Zener  tunnelling  time  is  derived  in  terms  of  the  applied 
unitom  electric  field  and  the  relevant  band  parameters.  The  analysis  shows  an 
eleelrio4ield-enhanced  broadening  of  the  excited  state  probability  amplitudes, 
thus  resulting  in  spatial  lattice  delocalization  and  the  onset  of  smearing  of 
disersM,  Stark  ladder  and  bancMo-band  transitions  due  to  the  presence  of  the 
electric  field,  in  addition,  the  velocities  of  a  Bloch  oscillation  will  be  observed  only 
for  the  electron  that  Is  initially  in  a  Bloch  state  before  Zener  tunnelling.  Further, 
the  influence  of  electnc  fields  on  resonant  tunnelling  structure  is  examined. 


1.  bt^oducllon 

Blod)  ekctroa  dynamics  in  a  homogeneous  electric  field 
has  been  a  sulg^  of  great  interest  dating  back  to  the 
earliest  applications  of  quantum  mechanics  to  solid  state 
physics  [1-3].  Even  more  recently,  as  modem  fabrication 
technologies  continue  to  drive  the  study  of  solid  state 
transport  into  the  nanometre  domain,  many  new  and 
interesting  questions  have  emerged  concerning  the  solid 
state  dynamics  and  quantum  transport  of  carriers  in 
‘band-engineered'  superlattices  and  tailored  periodic 
solids. 

In  this  paper,  a  novel  multiband  theory  of  Bloch 
electron  dyiuunics  in  homogeneous  electric  fields  of 
arbitrary  strength  is  presented.  In  this  formalism,  the 
dectric  field  is  described  through  the  use  of  the  vector 
potential;  in  this  regard,  this  work  builds  on  the 
methodology  previously  developed  by  one  of  the  authors 
(GJI)  and  co-workers  [1-3]  to  describe  solid  state 
dynamics  and  quantum  transport  for  Bloch  electrons 
in  an  apfritod  honu^eneous  dectric  field  of  arbitrary 
strength  and  time  dependence,  including  weak  scattering 
from  randomly  distributed  impurities  and  phonotu  [I], 
and  a  spati^y  localized,  inhomogeneous  electric  field 
[3].  In  addition,  multiband  coupling  it  treated  in  the 


Wigner-Weisskoirf’  (w-w)  [4]  approximation,  a  classic 
approach  for  describing  the  time  d^y  from  an  occupied 
quasistationary  state;  the  w-w  approxinution  allows  for 
the  analysis  of  the  long-time,  time-dependent  tunnelling 
characteristics  of  an  electron  transition  out  of  an  initially 
occupied  band  due  to  the  power  absorbed  by  the  electric 
field,  while  preserving  conservation  of  total  transition 
probability  over  the  comfdete  set  of  excited  bands. 

The  choice  of  the  vector  potential  gauge  leads  to  a 
natural  set  of  basis  functions  for  describing  Bloch  electron 
dynamics  in  a  homogeneous  electric  field.  A  basis  set  of 
localized,  elearic-field-dependent  Wannier  functions  and 
associated  envelope  equations  are  developed  to  accom¬ 
modate  the  inherent  localization  manifest  in  Bloch 
dynamics  due  to  a  relatively  strong  electric  field. 

Building  on  previous  methodology  [3],  use  is  made 
of  the  instantaneous  eigenstates  of  the  Hamiltonian 
describing  a  Bloch  electron  in  an  electric  field 

+  (1) 

where  I'(x)  is  crystal  potential,  Ao*  -  cE^t  is  the  vector 
•  potential  for  the  time-independent  homogeneous  electric 
field,  £o,  turned  on  at  initial  time  r  >  0.  The  solution 
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ABSTRACT 

The  F-X  scattering  rate  of  electrons  in  type-II  superlattices  by  optical-phonon  emission  is 
calfulated.  The  tight  binding  method  for  electronic  band  structure  and  the  dielectric  continuum 
modri  foe  phonons  are  used.  The  relative  strength  of  scattering  due  to  different  phonon  modes 
is  examined  fee  varying  si^crlattice  dimensions.  The  scattering  rate  is  highest  when  the  energy 
separation  between  the  r  and  X  levels  is  smallest,  and  decreases  qnlcUy  as  the  separatton  increases. 
It  is  found  that  the  strongest  scattering  ride  is  due  to  the  emission  of  AlAs  eowOwd  modes. 
Changing  of  parity  with  layer  thickness  and  Its  effect  on  scattering  are  discussed. 

1.  INTRODUCTION 

Recently,  there  has  been  considerable  interest  in  hot  carrier  dynamics  in  heterostmetures. 
One  impOTtant  issue  in  this  area  has  been  the  rdaxation  of  photoezdted  carriers  in  heterostme- 
tures  from  higher  energy  subbaads  to  the  lower  ones.  This  process  has  been  studied  extensivdy 
both  experimentally  and  theoretically  because  of  its  fundamental  physics  and  for  possible  device 
applicatiims.*  Most  of  the  work  up  to  now  has  been  on  type-I  heterostructures  whoe  the  bandg^ 
(done  material  is  entirely  nested  within  the  g^  of  another  material.  In  this  structure,  dectrons  and 
holes  are  both  localised  in  the  same  layer.  On  the  other  hand,  in  the  so-called  type-II  superlattices, 
the  hdes  are  localised  in  one  layer,  whereas  the  lowest  energy  dectrons  are  contained  in  the  other 
layer.^  The  wdl  studied  AlsGai.^As/GaAs  superlattices  can  be  made  type-II  by  appropriatdy 
choosing  the  layer  thicknesses  and  alloy  composition.  For  ezany^le,  the  GaAs/Al^  nqierlatticcs 
with  thick  AlAs  layers  are  known  to  be  type-II  for  GaAs  layer  thicknesses  less  than  or  equal  to 
35  A  (12  mondayers).^  Type-I  to  type-n  transitions  can  also  be  achieved  by  i^lication  of  external 
forces  such  as  an  electric  fidd  or  hydrostatic  pressure.  In  the  ease  of  type-II  structures,  the  lowest 
conduction  band  levd  is  not  in  the  GaAs  but  in  the  Al^Gai.^As  layer.  In  type-I  heterostmetures, 
dcetrtms  photoexdted  firom  the  valence  band  to  the  conduction  band  relax  to  lower  energy  levds 
within  the  same  layer.  However,  in  tjrpe-H  stractures,  the  dectrons  exdted  to  a  direct-gap  energy 
^d  (F)  in  COM  layer  then  can  relax  to  an  indirect-gap  energy  levd  (X)  in  the  a4}acent  layer. 
^  interesting  case  because  electron  relaxation  h^tpens  with  r^-space  transfis.  This 
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ABSTRACT 

Tli^  fTiwIltniiten  d— <!wlitiig  patutlal  tJ  eemSiamA  aeaoatic  phoMOna 

with  carxien  is  dariTod  by  qoantiang  tha  approptiata,  aq^athnantany-aaiffiad  appraadmate  cona- 
p— AimI  aganafelg  ptimuMi  ;ti  »  qwMitrww  mttm.  Tho  Scattering  ntO  doa  tO  tlw 

dafinmatioa  potential  interactkm  is  calculated  fiv  a  tango  of  qaantnm  wire  dfanensiana. 

1.  INTRODUCTION 

Idany  proposed  implications  of  maaoacepic  electronic  stmctnxaa  inaohw  carxiar  tranm<vt  at 
low  tenoperatnias  and  low  carxiar  anergias;  fraqoaatly,  the  la^ma  of  intarast  is  one  wbera  dhnan- 
tfamal  confinement  modifies  the  phase  space  snbatantially.  Li  this  low  tenmecatma,  low  energy 
repnie,^~*  aconstic  phoncns  plsqr  an  enhanced  nde  in  carxiar  scattering  and  may  dominate  over 
the  scattering  of  carriers  by  optical  phoncns.  Fortharxnore,  fat  nanoarale  stmetnres  it  is  poaalMe 
that  phase  tpatm  restrictions  may  weakm  or  fiarbid  imtical  phonon  scattering  procaaaea  that  would 
iMuniMiiy  Wi  Kiilfc  atwigttw-  In  lucont  ysacs,  them  hss  besB  on  extcBsiTe  Uterataxu  OB  tho 

role  of  dimensional  confinement  in  modifyiiig  longitiidiiial  optical  (LO)  phonon  modus  and  their 
interactions  with  charge  carxkrs  in  "■«*‘**^t*  and  masoaropic  semicondnctor  stmctnxes  (sea,  fior 
^“*^2  and  the  anmerons  pqm*  tfftrtnced  therein).  Howem,  thoe  are  xeiatie^ 
Urn  treatments  with  the  xok  of  dhnenaianal  confinement  in  modi^l^  aconstic  phonon 

modes  and  thdr  interactkns  with  charga  eaxxiers.*~^  WhQe  there  is  an  extendve  literature  on  the 
theory  of  aconstic  modus  in  conventional  wavegnidas,  resonators  and  xdated  stmctnres,  ftw  effiarts 
have  been  reported  to  fommiate  a  theory  of  aconteic  phonorn  in  nanoarale  stmctnres  adiece  both 
phonon  m«iii  ^  qnaatnm  treatment  of  ph******w  aorxnalimtion  are  essential} 

Constantinoa  has,  however,  dlscnstad  the  acoastie  phosion  modea  in  a  cylindrical  po¬ 

lar  sendcondncter  quantum  wire.^  The  need  fior  such  theoretical  treatments  has  been  underscored 
recently  by  eiperimental  studies  providing  both  direct  and  indirect*^  evidence  of  the  fanportaacs 
of  acoortic  phonon  confinement  in  reduced  «***"*"«^»t  dectronlc  structures. 
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The  r-JT  tnuuitioa  noe  for  eleetroaa  in  type-ll  luperlattket  it  calculated  for  the  case  of  optical* 
plMoon  wnwrion.  The  tiffat-brndint  method  for  ekctronic  band  ftructure  and  the  dwlectiic  continuum 
awdd  for  phonons  are  uaed.  The  relative  atrenfth  of  mattetinf  due  to  different  phonon  modes  is  eiam- 
ined  for  varying  saperiattice  dimensions.  The  scattering  rate  is  highest  when  the  energy  separation  be¬ 
tween  the  r  and  X  levels  is  smallrat,  and  decrewes  quickly  u  the  separation  increases.  It  is  found  that 
the  strongest  scattering  rate  is  due  to  the  emission  of  AJAs  confined  modes.  Changing  of  parity  with  lay¬ 
er  thicknem  and  its  eSbet  on  scattering  are  discussed. 


L  INTBODUenON 

Recently,  there  has  been  considerable  interest  in  hot- 
carrier  dynamics  in  heteroatructures.  One  important  is¬ 
sue  in  this  area  has  been  the  relaxation  of  photoexcited 
carriers  from  higher-energy  subbands  to  the  lower  ones. 
This  procem  has  been  studied  extensively  both  experi¬ 
mentally  and  theoretically  because  of  its  fundamental  na¬ 
ture  as  wdl  as  for  possiUe  device  applications.''^  Most 
of  the  work  up  to  now  has  been  on  type-I  heterostruc- 
tutes  where  the  band  gap  of  one  material  is  entirely  nest¬ 
ed  within  the  band  gi|>  tX  another  material.  Consequent¬ 
ly,  eketrons  and  bdes  are  both  localized  in  the  same  lay¬ 
er  in  these  structures.  On  the  other  hand,  in  the  so-called 
type-U  superiattices,  the  bedes  are  localized  in  one  layer, 
whereas  the  kmest-energy  dectrons  ate  contained  in  the 
opposite  layer.*  This  is  due  to  the  staggered  band  align¬ 
ment  in  these  heterostroctures.  The  well-studied 
Al,Ga,_j,As/OaAs  superiattices  can  be  made  type-II  by 
appropriately  choosing  the  layer  thicknesses  a^  alloy 
compositkms.  For  example,  GaAs/AlAs  superiattices 
with  thick  AlAs  layers  are  known  to  be  type-II  for  GaAs 
layer  thicknesses  Im  than  or  equal  to  3S  A  (12  mono¬ 
layers).*  Type-I  to  type-II  transitions  can  also  be 
achieved  by  application  cX  external  forces  such  as  an  elec¬ 
tric  field  or  hydrostatic  pressure.  In  the  case  of  type-II 
structures,  the  lowest  co^uction-band  level  is  not  in  the 
OaAs  but  in  the  Al^Gsi-^^As  layer.  In  type-I  hetero¬ 
structures,  electrons  photoexcited  from  the  valence  band 
to  the  cc^uction  band  relax  to  lower-energy  levels 
within  the  same  layer.  However,  in  type-II  structures, 
the  dectrtms  excited  to  a  direct-gap  energy  level  (D  in 
oat  layer  then  can  relax  to  an  indirect-gap  energy  level 
(AT)  in  the  adjacent  layer.  This  is  a  very  interesting  case 
because  electron  relaxation  occurs  in  conjunction  with 
real-q»ace  transfer.  This  process,  which  is  normally  for¬ 
bidden,  is  now  possiUe  due  to  the  mixing  of  F  and  X 
states  by  the  snperlattice  potential  and  relaxation  of 
momentum  conservation  due  to  interface  disorder  Qateral 
mixing).  Several  experimental  resultt  related  to  F-JT 
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transfer  are  available  in  the  literature.*'^  In  spite  of  the 
spatial  charge  transfer,  the  relaxation  rate  can  be  very 
high.  Ejqierimental  evidence  has  been  presented  that 
F-Jf  electoon  transfer  occurs  by  the  emission  of  long- 
wavelength  optical  phonmu.*  However,  a  detailed 
theoretical  analysis  for  the  T-X  transfer  in  type-H  struc¬ 
tures  has  not  yet  been  given. 

In  this  p^>er,  are  calculate  the  F-Jt  scattering  rates  due 
to  optical-phonon  emission  using  a  realistic  band- 
structure  model  and  show  that  this  is  an  important  mech¬ 
anism  for  the  F'jlf  relaxaticm  process  in  ty^n  superiat- 
tices.  The  paper  is  organized  as  follows.  In  Sec.  II,  we 
present  the  ti^t-binding  model  and  the  dielectric  contin¬ 
uum  model  used  in  the  calculation.  Evaluation  of  matrix 
elements  is  explained.  In  Sec.  Ill,  we  present  and  discuss 
our  results  of  electnm  relaxation  rates.  Rnally,  in  Sec. 
rv,  a  summary  ir  given. 

n.  FORMULATION 

It  is  well  known  that  the  electron-optical-phonon  in¬ 
teraction  in  low-dimensional  systems  can  be  altered 
strongly  due  to  the  confinement  of  the  carriers  and  the 
confinement  of  the  phonons.  Usually,  a  single-band, 
spherical  effective-mass  model  is  used  for  the  description 
of  the  confined  carrier  states.  It  is  commonly  assumed 
that  each  of  these  states  is  derived  from  bulk  states  of  a 
given  symmetry  only  (e.g.,  F,  X,  or  L),  and  leveb  derived 
from  different  bulk  states  do  not  interact  with  each  other. 
In  a  superlattice,  momentum  in  the  growth  direction  is 
not  conserved  due  to  the  discontinuities  in  the  snperlat¬ 
tice  potential  in  this  direction.  Therribre,  bulk  states  of 
the  same  energy  but  different  symmetry  can  couple  with 
each  other.  For  a  superlattice  grown  in  the  z  direction, 
the  sixfold  degeneracy  of  the  X  valley  is  removed,  with 
the  formation  of  an  X^  doublet  having  momentum  along 
the  growth  axis,  and  an  Xg^y  quadruplet,  with  momentum 
in  the  plane.  The  and  L  minima  are  not  coupled  to 
the  F  valley  due  to  the  conservation  of  lateral  momen¬ 
tum.  Therefore,  only  the  F  and  Xg  valleys  are  coupled 
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EleetnM-aoointiciilioiKNi  icaftffine  in  a  reettngular  quantum  wire  ia  atudkd.  The  Hamiltonian 
rtwrriWnj  the  defocmatkm-poiaatial  interaction  of  conSned  acouatk  phonona  with  carrieta  ia  derived  by 
qnantiriag  die  appropriate,  eiperimentaUy  veriied  approximate  romprfiaional  acouatic-phonon  modea 
ia  a  free-etandin(  reetansular  quantum  wire.  The  acattering  rate  due  to  the  defcnmation-potential  in- 
tcractkm  ia  obtained  for  OaAa  quantum  wirea  with  a  range  of  croaa  aectional  dimenaiona.  The  reaulta 
demonatrate  that  a  proper  treatment  of  confined  acouatic  phonona  may  be  fiaential  to  correctiy  model 
electron  acattering  ratea  at  low  eoergiea  in  nanoacale  atructuiea. 


L  INIROOUCnON 

A  number  of  piopoaed  appUcationa  of  meaoacopic  elec¬ 
tronic  atmettirea  invtdve  carrier  transport  at  low  temper- 
atnrea  and  low  carrier  energiet;  frequently,  the  regime  of 
interwt  ia  one  where  dhnenakmal  confinement  modifies 
the  phase  apace  substantially.  It  ia  weU  known  that  in 
tUa  low-temperature,  low-energy  regiine,*'^  acoustic 
l^ionana  play  an  enhanoed  role  in  carrier  acattering  and 
may  dominate  over  the  acattering  ci  carriers  by  optical 
phonons.  In  addition,  in  nanoacale  structures  it  is  possi¬ 
ble  that  friiaseepace  lestrictiona  may  weaken  or  forbid 
optical-phonon  acatteiiag  processes  that  would  normally 
dominate  in  bulk  atmeturea.  Recently,  there  has  been  an 
extensive  hteratnre  on  the  role  of  dimensitmal 
confinement  in  modifying  kmgttodinal-optical  (LO)  pho¬ 
non  modes  and  thdr  interactions  with  charge  carriers  in 
nanoacale  and  mesoscopic  semiconductor  structures  (see, 
for  exanqile.  Reft.  7-12  and  the  numerous  pliers  refer¬ 
enced  therein).  On  the  other  hand,  there  are  relatively 
few  treatments  dealing  with  the  role  of  dimensional 
confinement  in  modifying  acoustic-phonmi  modes  and 
their  interactions  with  charge  carriers.^~^'^’'*  In  spite  of 
the  fact  that  there  is  an  extensive  literature  on  the  theory 
of  acoustic  modes  in  conventional  waveguides,  resona¬ 
tors,  and  related  structures,  few  efforts  have  been  report¬ 
ed  on  formulating  a  theory  of  acoustic  phonons  in  nano¬ 
acale  atnaetnrea,  ediete  both  phonmi  confinement  and  a 
quantum-mechanical  treatment  ot  phonon  normalization 
are  essential.  The  necessity  for  such  theoretical  treat¬ 
ments  has  been  demonstrated  recently  by  experimental 
studies  |»oviding  both  direct  and  indhect^'^  evidence  of 
the  importance  of  aconstic-phmion  ermfinement  in  re¬ 
duced  dimensional  electronic  structures. 

In  dus  pqper,  we  have  obtained  the  normalized  expres¬ 
sions  fbr  acoustic  phonons  confined  in  a  free-standing 
rectangular  quantum  wire  by  aiqiropriately  quantizing 

SO 


the  acoustic-phonon  displacements.  As  is  well  known, 
there  are  no  exact  solutions  for  the  complete  set  of  pho¬ 
non  modes  for  a  rectangular  quantum  wire;  nevertheleas, 
as  for  the  case  of  LO  phonon  modes,*  the  approximate 
modea  |»esented  in  this  woric  provide  simple  and  naeftal 
expressions,  which  are  well  suited  (or  modeling  the  in¬ 
teraction  of  carriers  with  acoustic  phonona.  As  a  basis 
for  investigating  the  role  reduced  dimensionality  on 
the  coupling  between  acoustic  phonons  and  carriers,  we 
have  formulated  the  interaction  Hamiltonian  for  the  de¬ 
formation  potential  associated  with  confined  acoustic- 
phonon  modes  in  rectangular  quantum  wires.  The  result¬ 
ing  scattering  rates  (baaed  on  the  gdden  rule  ^iproxima- 
tkm)  are  compared  with  those  obtained  fnmi  Ae  bulk- 
phonon  modes.  For  numerical  calculations,  OaAs  is  used 
as  the  material  of  choice  throughout  this  study. 

a  QUANTIZATION  OF  COMPRESSIONAL 
ACOUSTIC-PHONON  MODES 
FOR  A  RECTANGULAR  QUANTUM  WIRE 

The  compressional,  or  dilatational,  acoustic-phemon 
modes  in  free-standing  rods  of  rectangular  cross  section 
have  been  examined  both  experimentally'*  and  theoreti¬ 
cally'**'^  by  Morse  in  an  extended  study.  Morse  has  de¬ 
rived  an  approximate  set  of  hybrid  compresskmal,  or  di- 
lataticmal,  acoustk-phtmmi  modes.'^'^  which  are  found 
to  accurately  approximate  the  experimentally  observed 
modes  over  a  iride  range  of  conations.'*  ^lecifically, 
Morse  has  found  that  the  approximate  h]rbrid  modes  d^ 
rived  by  ■—«««"«";  sqiarable  boundary  conditums'**'^ 
have  simple  analytical  representations  1^  provide  con¬ 
venient  approximations  for  the  rectangulu'  geometry 
when  the  cross-sectional  dimenakms  have  aspect  ratios 
approximately  2  or  greater.  For  smaller  aqiect  ratios 
(i.e.,  close  to  I),  Morse  has  argued  correctly  that  it  is 
necessary  to  turn  to  numerical  solutimis  since  exact 
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